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Abstract 

Several basic problems of the theory of quantum phase transitions are reviewed. 
The effect of the quantum correlations on the phase transition properties is con- 
sidered with the help of basic models of statistical physics. The effect of quenched 
disorder on the quantum phase transitions is also discussed. The review is performed 
within the framework of the thermodynamic scaling theory and by the most general 
methods of statistical physics for the treatment of phase transitions: general length- 
scale arguments, exact solutions, mean field approximation, Hubbard-Stratonovich 
transformation, Feynman path integral approach, and renormalization group in the 
field theoretical variant. Some new ideas and results are presented. Outstanding 
theoretical problems are mentioned. 

1. Introduction 

1.1. Historical notes 

The experimental and theoretical research of quantum phase transitions and quantum 
critical phenomena is a branch of the statistical physics of a rapidly growing importance for 
the explanation of essential features of low dimensional fermion and spin systems, dilute 
Bose fluids, superconductors, quantum Hall systems, ferroelectrics (HEHHHUIEIIHIEE]- 

The quantum statistics has a substantial influence on the critical behaviour near various 
(multi) critical points of low temperature and zero temperature continuous phase tran- 
sitions as well as on thermodynamic and correlation properties near equilibrium points 
of low-tempetarure first order phase transitions [HJ E3 HTj . This involves the quantum 
statistical physics in the field of a particular type of phase transitions - the quantum 
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phase transitions. The thermodynamics and the correlation phenomena near such phase 
transitions are essentially affected by quantum mechanical effects. 

It seems convenient to begin this review with brief historical notes. Although the inter- 
est in quantum effects on the critical behaviour in many-body systems dated from the 
dawn of the quantum statistical physics, the theory of quantum phase transitions and, 
in particular, the theory of continuous quantum phase transitions (alias, quantum criti- 
cal phenomema) received a substantial development after the late sixties of the previous 
century owing to the application of fruitful ideas of universality and scaling by renormal- 
ization group (RG) methods; see, e.g., Refs. [TJ H21 EH EH E3 EH] • The contemporary 
theory of these phenomena has two distinct periods of development. 

The first, classical period, began in the seventies of the previous century with the pio- 
neering papers by Pfeuty and Elliott [E| , Rechester [THj , Young [TH] , and the particularly 
important work of Hertz (SHj. This period of a relatively quiet research based on quantum 
field and statistical methods, including RG, continued till 1986-1987 when the number of 
papers on quantum critical phenomena has abruptly increased due to the great interest in 
heavy- fermion [Hj and high-temperature superconductors [21], low- dimensional magnetic 
systems jO], quantum Hall effect metal-insulator transition problems 01221, an d dilute 
Bose gases and Bose fluids P El 123 123 ■ 

The new period of research [21 E3 HI El EH IE] after 1986-1987 seems to be much more 
intensive but one should not forget that this renewed research in an extended area of 
problems and their applications relies very much on the classical results obtained during 
the first period. We shall partly consider some classical results together with related 
developments and applications accomplished up to now. 

1.2. Scope and aims of the review 

Numerous real systems, for example, magnets, ferroelectrics, quantum Hall systems, 
dilute Bose gases, low-dimensional Josephson-j unction arrays, heavy-fermion and high- 
temperature superconductors exhibit phase transition lines which extend to low and up 
to zero-temperature phase transition points, namely, to temperatures of a strong quan- 
tum statistical degeneration P0E1II1E1E1II1E1I2I1I221I2II1I2I]- In this temperature 
region the quantum statistical correlations can, in certain cases and under certain exper- 
imental conditions, produce observable effects on the low- and zero-temperature phase 
transition. The task of the theory is to predict and describe the quantum effect on the 
thermodynamic and correlation properties near such phase transitions. This difficult task 
remains among the outstanding problems of the quantum statistical physics, despite the 
remarkable success achieved in this field of research during the last 30-35 years. 

We shall not expand our discussion over the great variety of low-temperature phase tran- 
sitions. Rather we shall focus our attention on several aspects of continuous and, mainly, 
second order quantum phase transitions, i.e. quantum critical phenomena. For our aims 
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we shall use general thermodynamic arguments and esults for basic and relatively simple 
statistical models. 

Despite of the fact that in real systems at low temperatures the first order phase transitions 
occur more frequently, quantum phase transitions of first order are less investigated and 
few results are available for them. Perhaps, this is due to the fact that the ordering 
phenomena and thermal fluctuations at equilibrium points of discontinuous (first-order) 
phase transitions have finite length scales and the scaling methods cannot be applied. 
That is why the quantum fluctuations at first order phase transitions should have a 
stronger effect on the phase transition properties than at points of continuous phase 
transitions. This interesting topic is also beyond the scope of the paper. 

Our attention will be focussed on the quantum behaviour of Bose gases which reveals 
the main features of the quantum phase transitions in various systems. Furthermore, our 
consideration is closely related to the problem of Bose-Einstein condensation (BEC) of 
noninteracting bosons (ideal Bose gas, shortly, IBG) and the superfluidity of interacting 
bosons (nonideal Bose gas; NBG) EH HI]. We shall be particularly interested in the 
properties of usual and disordered Bose fluids in a close vicinity of the A-point. General 
problems of quantum phase transitions are also reviewed. Along with this review we shall 
present several new aspects of the quantum phase transitions. The new ideas and results 
will be summarized in our concluding remarks. 

In Sec. 2 we perform a phenomenological investigation of several topics concerning the 
change (crossover) of the usual classical ("high-temperature") critical behaviour when 
the critical temperature is lowered to the range of the quantum degeneration up to 
zero. The lowering of the critical temperature is referred to as "low-temperature" and 
"zero-temperature" limiting cases of the critical behaviour. Numerous experimental and 
theoretical studies indicate that a rather nontrivial and general "high-low" temperature 
crossover (HLTC) exists, namely, that the high-temperature and low-temperature crit- 
ical phenomena are quite different from each other. The first fundamental question is 
whether this HLTC is a result of low-temperature and zero-temperature limits themselves 
or is produced by quantum effects. 

The problem can be solved provided both the universal and nonuniversal properties of 
the low-temperature critical behaviour are thoroughly investigated. Note, that some par- 
ticular properties of the system or special experimental conditions can quell the quantum 
correlations and then the low-temperature critical behaviour will remain totally or par- 
tially classical up to the zero temperature. Therefore, the quantum critical behaviour is a 
particular case of low-temperature and zero-temperature critical behaviour. The classical- 
to- quantum dimensional crossover (CQC) [T7[ 01201 is also discussed as a particular case 
of HLTC. For our aims in Sec. 2 we apply a new scaling scheme for quantum systems |5J 
which is an alternative of the standard quantum scaling, developed for the first time in 
Ref. (see, also, Refs. |01 El I2U) 
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In Sec. 3 IBG at constant density and constant pressure is considered. The general 
methods for the investigation of systems of interacting bosons are reviewed in Sec. 4. In 
Sec. 5 the main RG results for NBG are discussed. In Sec. 6 the transverse Ising model 
(TIM) is considered with the help of the mean-field (MF) approximation, the consideration 
of the Ginzburg critical region PQE], and RG. The results in Sees. 3, 5 and 6 confirm the 
scaling description outlined in Sec. 2. Disorder effects are discussed in Sec. 7. In Sec. 8, 
the notion about the break down of universality for quantum critical phenomena and a 
classification of these phenomena based on this notion are established. The concluding 
remarks are presented in Sec. 9. We presume that the reader is well acquainted with the 
basis of the scaling theory of phase transitions and RG [TJ [T3J EH UJJ U3] • 

There is an overwhelming amount of papers published in this field and we cannot enu- 
merate and discuss all of them within a manageable length of this review. We have tried 
to present the most important original papers, comprehensive review articles and books. 



2. Phenomenology and quantum phase transitions 

2.1. General criterion for quantum effects on phase transition 

The quantum effects (alias, quantum correlations, often called also quantum fluctuations 
and even quantum pseudo- interactions) due to the overlap of particles wave functions exert 
an influence on the thermodynamic and correlation properties near the phase transition 
points in many body systems provided the de Broglie thermal wavelength 

^ ■ (SO'- 

is greater than the correlation length 

i{T) = UT c )\t(T)r (2) 

of the thermal fluctuations: 

q = | > 1. (3) 

In Eqs. (1) - (2), 

0<|t(T)| = ^^«l, (4) 

defines a broad vicinity of the (multi) critical point T c , in which the phase transition 
phenomena occur, £ = £(T C ) is the so-called "zero temperature correlation length," the 
critical exponent v > describes the behaviour of the correlation length £ (T) in the phase 
transiton region (4), 9 > is the exponent for the thermal wavelength. We have to note 
that usually, 9 = 1/2 [TU] but in order to comprise all possible quantum statistical models 
we consider 9 > 0; see, e.g., [TJ. In Eq. (1), the parameter m denotes either the real 
particles mass or an effective mass m e g ~ h 2 /c(g) of composite bosons (boson excitations) 



4 



which represents the effect of some interaction constant g; below the suffix "eff" will be 
omitted (see also Sees. 3.2 and 4.1). 

We must keep in mind that the correlation length £ (T) describes only classical phenonena. 
The length £ is a nonuniversal quantity and can be specified only for concrete systems 
(models) which means that it cannot be presented in a concrete mathematical form before 
the choice of the concrete model (Hamiltonian) is done and, in particular for 
example IBG, after the thermodynamic analysis of the system is fulfilled. That is why, 
at this general stage of consideration we do not give a mathematical formula for £o- But 
we must emphasize that £ depends on intrinsic parameters of the system and, especially, 
on the critical temperature T c . The term "zero-temperature correlation length," that is 
sometimes used for the scaling amplitude £ of £(T) is an indication of the simple fact 
that t(0) = 1 and, hence, £(£ = 1) = £o- 

It will be shown in Sec. 4.3 that 9 = 1/z, where z is the so-called dynamical critical 
exponent which describes the intrinsic (quantum) dynamics of the system [20] ■ It can be 
stated as a theorem that the relation 6 — 1/z is valid always when the dynamics of the 
quantum system [20] is not influenced by other time- dependent phenomena [T2"l I25j. 

It has been already mentioned that T c is a (multi) critical temperature, i.e., it denotes 
the phase transition point of quantum or classical phase transitions of second and higher 
order, i.e., continuous phase transitions. For the first order phase transitions, however, 
Eqs. (2) and (4) can be also used. In this case T c is a characteristic temperature of the 
system near the equilibrium phase transition temperature T eq ^ T c , (\T — T eq \/T eq ) = 
t eq (T) ~ \t\ -C 1. The phase transition phenomena, including the metastable states are 
located in the temperarute domain \t eq (T) ~ t eq (T c )\ pp. The infinitesimally small vicinity 
[(t(T) — > 0] of critical points T c of continuous phase transitions is often referred to as an 
"asymptotic critical region" and the behaviour in this small domain is usually called an 
"asymptotic critical behaviour." 

Phase transitions which satisfy the criterion (3) in the vicinity of their equilibrium phase 
transition points T c are called quantum phase transitions. This criterion has been intro- 
duced for the first time by M. Suzuki for quantum critical phenomena (phenom- 
ena at continuous quantum phase transitions) and here we find reasonable to extend 
it to all quantum phase transitions. The condition (3) is very similar to the respec- 
tive one for quantum fluids, for them there is a requirement for the thermal length 
A to exceed the mean interparticle distance (the lattice constant in crystal bodies): 
a = p x l d = (N/Vfl d OH N is the particle number, V = (Li...L d ) is the vol- 
ume of the system, and d > is the spatial dimensionality. 

The criterion (3) is a direct result of the fundamental notion that the phase transition is 
a (quasi)macroscopic phenomenon and, hence, phenomena at length scales shorter than 
the length scale £ of the classical fluctuations are irrelevant [26J. It is confirmed by the 
available studies of quantum statistical models. The criterion (3) yields both the necessary 
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and the sufficient conditions for the appearance of a quantum phase transition but the 
circumstances, under which this can actually happen remain hidden in its generality. 

It seems at first glance that the condition (3) can be more easily fulfilled for first order 
transitions, where £(T e9 ) < oo. However, the answer of the question whether the in- 
equality (3) can be satisfied even in this case depends on the values of the parameters £ 
and Ao = (2ith 2 jmhsf = XT 9 for the concrete substance of interest. The same problem 
exists for the continuous phase transitions. Hereafter we shall discuss continuous phase 
transitions, where £(T) is always infinite. 

2.2. Quantum critical region 

The vicinity (4) of the critical point T c , where the phase transition phenomena occur and 
the Landau expansion of Gibbs free energy in powers of the fluctuation order parameter 
(f)(f) is valid, includes both the MF domain of description and the Ginzburg critical region 
of strong fluctuations very close to T c [TJ[TT]. This vicinity of T c , £(T) 3> £q(T c ), where the 
classical (thermal) fluctuations and ordering phenomena occur, hereforth will be called 
"classical transition region" or, shortly, "transition region". For some types of phase 
transitions the transition region may be defined by \t(T)\ < 1 rather than by the strong 
inequality (4). 

While the correlation length £(T) does not describe the MF domain of ordering and, 
moreover, this quantity cannot be defined within the standard MF approximation, the 
inequality £(T) 3> £,o(T), equivalent to t(T) <C 1, does include the temperature domain of 
the phase transition phenomena and can be used in our discussion. We shall distinguish 
between the Ginzburg critical region |IJ and the much larger classical transition region 
defined above. The subdomain of the transition region, where the quantum fluctuations 
have an essential influence on the phase transition properties will be called a "quantum 
transition (sub) domain." At a next stage of the present consideration, this picture will be 
generalised to include additional thermodynamic parameters. 

From Eq. (1), the condition (3) and £o(T c ) > ao - a condition which is always satisfied in 
real systems, we see that a quantum critical phenomenon may occur only in the regime of 
quantum degeneration of the system (A > a), i.e., near low- and zero-temperature critical 
points, where conditions (3) and (4) are simultaneously satisfied. Here we accept as a 
definition that we shall use the term "quantum critical phenomenon" to indicate the fact 
that the quantum effects do penetrate in the transition region |t| <C 1. 

Quantum critical phenomena, that is, quantum effects in the "transition region", may 
occur in condensed matter systems and gases, where the critical line T C (X) extends over 
the low-temperature region to the zero-temperature critical point T C (X ) = by variations 
of an additional thermodynamic parameter X. Depending on the particular system, the 
parameter X may represent some auxiliary intrinsic interactions which suppress the main 
interaction responsible for the phase transition, or other physical quantities (density, 
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Figure 1: (a) Low temperature part of a critical line with zero temperature critical point 
T C (X Q ) = 0. The shaded part (a — X — b) of the transition region, marked by the lines 1 
and 2, corresponds to a low temperature classical behaviour, (b) Low temperature critical 
line with X = 0; domains (a — X — b) and 1 — — 2 coincide. 

pressure, concentration of impurities, etc.) having the same effect; see Fig. 1, where such 
critical lines are depicted for Xo > and Xo = 0. The slope of the curve T C (X) and the 
fact that it is tilted to the left are irrelevant to the present discussion. In certain systems, 
the critical line T C (X) is tilted to the right and then the maximal value of T c will be at 
some X > X (see Sec. 6.2). 

We should keep in mind that Eqs. (2) and (4) describe only the temperature dependence of 
the characteristic length £. In general, the correlation length £ depends on the parameter 
X as well: £ = £(T, X). The dependence of thermodynamic and correlation quantities on 
X must be considered on the same footing as the dependence on T. The reason is that, 
as seen from Fig. 1, the phase transitions at T C (X) or, equivalently, at X C {T), can occur 
both by variations of T at fixed X [T -driven transitions) and by variations of X at fixed 
T (X -driven transition). Phase transitions produced by both T- and X-variations are 
also possible. We suppose that their main properties can be clarified by the investigation 
of T- and X-transitions. 

These notes are important for the evaluation of that part of the transition region where 
the quantum effects penetrate, i.e., in the quantum transition region. The inequality (4) 
defines the transition region of the critical phenomena produced by T-transitions and, 
hence, the temperature transition width [T>(X) — T<(X)] along the T-axis satisfies the 
condition [T>(X) - T<(X)] < T C (X) for each X = X c ; here, T>(X) > T<(X). As 
T C (X) — > for X — > X , the temperature transition width vanishes at the point (T, X) = 
(0, Xo). The shape of this temperature transition region is similar to the form of the 
shaded domains a - X - b and 1-0-2 in Figs. la,b. 
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However, the transition region along the X-axis is produced by X-transitions and this 
effect should be considered as well. That is why the notion for the transition region 
should be generalized to include the X— driven phenomena, too. At any given T = T c , 
£(T C , X) — > oo for X — > X c because of the general criticality assumption: £(T C ,X C ) = oo. 
Therefore, one can represent £(T C , X) in a general scaling form, namely, £(T C , X) = £ot~ u , 
where t(X) — \X — X c |/X c is the X— distance from X c , v > is the respective critical 
exponent, and £ = £o(X c ) = £[T C ,£(0)] > is the scaling amplitude. Now one can 
define the X-transition interval [X>(T) - X<(T)] < X C (T) for each T = T c ; here 
X>(T) > X < (T). This interval will vanish, only if X C (T) = at some critical temperature 
(T = T c ). The latter case is illustrated in Fig. la by the high-temperature critical point 
(T c ,0), and in Fig. lb, where X c (0) = 0. 

The total transition region (T, X)t in the vicinity of the critical line T C (X) includes both 
T— transition and X— transition widths and can be represented as a two-dimensional 
domain of shape similar to the form of the domain confined between the lines 1 and 2 
in Figs. la,b. The picture of the transition region (T,X) t outlined above will be valid 
provided the scaling amplitude Co(T c ) is finite for all T c > 0. Remember that (T, X) t will 
be a classical transition region till the quantum effects are ignored in our consideration. 

The next question is whether the quantum correlations, if properly taken into account, 
will affect the thermodynamic properties in the transition region (T,X) t . When the 
transition region (T, X) t vanishes, as this is the case for some zero temperature critical 
points, the quantum fluctuations will fill up the whole nearest vicinity of the critical point 
and then we have an entirely quantum phase transition. The problem is nontrivial when 
the classical region (T,X) t exists. 

There are three scenarios: 

(i) The quantum effects are relevant in the whole region (T, X) t . Then the latter is a 
quantum transition region, where quantum phase transition phenomena occur. 

(ii) The quantum effects are relevant in a part (quantum subdomain) of the total transition 
region. This will be enough to include the phase transition in the quantum ones. 

(iii) The quantum effects are irrelevant anywhere in (T, X) t and, hence, the phase tran- 
sition is completely classical. 

The variants (i)-(iii) are shown in Fig. la,b. The shaded domain a - Xq - b in Fig. la 
represents the classical subdomain of the total transition region (T, X) t , where the classical 
(thermal) fluctuations dominate and the phase transition is completely classical. The 
unshaded domain between the lines 1 and 2 is the quantum transition subdomain. Here 
one can find both cases (i)-(ii) depending on whether T > or T = and, whether one is 
interested in T— or X— driven phase transitions. In Fig. lb the classical domain a - X - 
b coincides with the total transition domain 1-0-2 which is an illustration of the classical 
variant (i). 
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In Sees. 2.3-2.4 we shall investigate the variants (i)-(iii) with the help of general scaling 
arguments. For a methodical convenience, we shall discuss the T— transition region along 
the T-axis and the X-transition region along the X-axis separately, although they are 
a result of one and the same reason - thermal and quantum fluctuation and ordering 
phenomena. 

2.3. T— driven transitions 

According to Eqs. (1) and (2), £ — > oo for T — > T c > 0, but A remains finite. Therefore, 
the close vicinity of finite temperature critical points T c > will always exhibit a classi- 
cal behaviour. In particular, this is true for the so-called asymptotic critical behaviour 
corresponding to the infinitesimally narrow distance (|t| — > 0) from T c . 

Using the criterion (3) as well as Eqs. (1) and (2), it is easy to show that quantum critical 
phenomena will occur, i.e., the quantum effects will penetrate in the (temperature) critical 
region \t(T)\ < 1, if 

T c {io/\) l/u < \{T-T C )\ « T c . (5) 

The conditions (5) are well defined for T c > 0, including the zero-temperature limit 
T c — > 0, provided T — > too, so that \t\ < 1 is satisfied. The size of the quantum 
subdomain (AT)q = \T* — T C \ above and below T c is given by (AT)q = Tc 1+e ^ u \^ / ^o) 1 ^ , 
where A = \(T)T 9 . The size of this temperature interval depends on the behaviour of 
the function £,o(T c ). 

Obviously, the quantum portion (5) of the transition region will exist, if 

£o(T c ) < A(T), (6) 

which corresponds to a moderate (A > a) or strong (A ^> a) quantum degeneration 
for (£ > a) and (£ 3> a), respectively. The condition (6) can be written in the form 
T° < [Xo/UT c )}. 

The inequality (6) may bring more information about the quantum criticality correspond- 
ing to T-driven transitions, if we remember that within the present description (\t\ <C 1 or 
\t\ < 1), we can freely substitute T with T c in all formulae except that for t(T). Conversely, 
the substitution of T c with T will be also allowed in the whole transition domain (4), if 
the nominator of t(T) is not affected. Besides, the latter variant of the theory seems to be 
more closely connected with the original form of the free energy PI , which is derived 
by standard statistical methods from the microscopic Hamiltonian of the system (see also 
Sec. 4). 

The further consideration depends on the way, in which we shall treat the nonuniversal 
length Co{T c ) > a. When T c is decreased to some extent, depending on the specific 
properties of a given system, £o(^c) grows and this is described by the relation 

UTc) = tooT-»° , (7) 
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where a < £00 < 00 is a new (low-temperature) scaling amplitude and < is < 00 is a new 
(low-temperature) exponent; z/ must be nonnegative number because of the criticality at 



The value of exponent vq and the range of temperatures, where the low-temperature de- 
pendence (7) will be valid, can be evaluated from the properties of system. The quantities 
£00 > o and uq > appear as a result of a gradual low-temperature crossover in the critical 
behaviour owing to T c lowering. 

The temperature, T, is not sufficient for the description of the low-temperature crossover 
and the entire investigation can be performed in terms of both relevant thermodynamic 
parameters T and X. The low temperature critical behaviour limit should be given 
by the ratio [T c (X)/T c (0)] <C 1 or, equivalently, by the ratio [g(X)/g ] <C 1, where 
9o = fi'(O) is the interaction responsible for the T-driven transition and g{X) describes the 
effective decrease of this interaction owing to the auxiliary parameter X; in certain cases, 
g(X) ~ T C (X). These remarks, involving the term "interaction" in our phenomenological 
analysis, are not restricted to systems of interacting particles only; see Sec. 3, where 
the phase transition is produced rather by the Bose statistics and global thermodynamic 
constraints than by direct interparticle interactions. 



Having in mind these notes we continue our analysis of the condition (6) by substituting 
A(T) with A(T C ) and (,o(T c ) with the scaling form (7): 



When the critical temperature is lowered enough, the condition (8) will be broken unless 



If the inequality (9) is satisfied, the quantum critical phenomena will occur at low- 
temperatures in certain part (5) of the temperature interval (4) and, moreover, for T c — > 0, 
they will prevail in the asymptotic critical behaviour [t(T) — > 0], too. 

If the criterion(9) is fulfilled, the classical fluctuations will be completely irrelevant to 
the critical behaviour at T c = 0. Their effect at T c ~ (extremely low-temperature crit- 
ical points) will be restricted in a negligibly narrow vicinity of T c , which is practically 
inaccessible to experiments. In both cases, T c ~ and T c = 0, the experiment should 
observe quantum critical phenomena only. For — vq, the quantum correlations will have 
an effect on the pre-asymptotic low-temperature critical behaviour outside the small dis- 
tance T c (£qo/\o) 1 / u from T c , provided (£00/^0) < 1> i- e -> m case °f convenient nonuniversal 
properties of the particular system. Quantum critical phenomena will not exist at all, if 
6 < vq. In this case the critical behaviour ramains classical. 

The same results can be obtained by an alternative consideration which is consistent with 
our discussion about the possible change of factors T with T c , and vice versa. The scaling 



low T c . 




(8) 



9 > Vq . 



(9) 
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law (2) has the following alternative definition: 

eco = uT)\t(T)r, (io) 

where t(T) = (T — T c )/T. The quantity £o(T) defined by Eq. (10) can be called "zero— T c 
correlation length"; i(T c = 0) = 1. In the high-temperature and, to some extent, in the low 
temperature range of critical temperatures T c , Eq. (10) will give the same leading scaling 
dependence as Eq. (2) because the corrections to the scaling are of the order 0(\t\ 1 ^ u ). 
In the limit T — > however the quantity t(T) is not singular, whereas the correlation 
length £(T) should tend to infinity at T c = 0. This means that the decreasing of the 
temperature will unavoidably yield a HLTC. The latter consists in a gradual temperature- 
driven transformation of the t(T)— singularity in Eq. (10) to a new singularity - the 
divergence of the amplitude £o(T) at T c = 0. 

Let us suppose that £o(T) obeys a scaling law with respect to T having the form (7). Then 
one can immediately rederive a condition for temperature T which will be identical to 
the condition (8) for T c and, hence, the criterion (9) will be straightforwardly confirmed. 
Note, that the new (low temperature) singulatiry may arise with respect to the vasiable 
t(X) rather than towards T (see an example in Sec. 6.2). 

The singularity (7) is a mere indication for an existence of a scaling law for £(T) of the 
type £(T) ~ £o{T) = £oo/T uo corresponding to zero critical temperatures (T c = 0). In 
the low-temperature limiting case, the critical behaviour singularities are developed as 
singularities of the scaling amplitude while at high temperatures this scaling amplitude is 
a slow varying function of temperature. It is worth noting that the present phenomeno- 
logical consideration, including the ansatz (7) for £,o(T c ) and the related one for £o(T), is 
supported by results from RG calculations for particular systems [28J and experiments on 
ferroelectric phase transitions (2H1 EDj • 

2.4- X -driven transitions 

In general, the properties of X— driven transitions are quite different from those of 
T— driven transitions. In particular, essential differences may be expected in the high- 
temperature (X ~ 0) and low-temperature (X ~ Xq) ranges of temperatures. By sub- 
stituting T with X and, of course, T c with X c in (2) and (4), one can perform the phe- 
nomenological analysis (Sec. 2.3) in terms of the variable X. This analysis directly yields 
a form of scaling law for the correlation length £ (T, X) which is identical to that related 
with T-transitions in the high-temperature region where X ~ 0. The only difference may 
come from the values of critical exponents which describe the scaling laws with respect to 
t(T) and t(X). There is no reason to suppose that the exponents towards t(T) and t(X) 
should be equal. 

In fact, the critical behaviour in the low-temperature limit described in Sec. 2.3 has a 
similarity with high-temperature X-transitions because X c ~ at high temperatures. A 
conformity of type T <-> X concerning the form of scaling laws and the way, in which 
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the low-temperature limit (for T-transitions) and the high-temperature limit (for X- 
transitions) are developed from the corresponding scaling amplitudes certainly exists. 
However, except for particular systems, a total correspondence including the values of 
respective critical exponents and scaling amplitudes cannot be expected; note, that the 
length A(T) has no analog in terms of X. In contrast, there is no such similarity between 
the low-temperature T- and X-transitions. 

For X-transitions, the criterion (3) gives 



\t(X)\ > 



l/v 

A(T) 



(11) 



where all quantities are defined by the change of T with X in (2) and (4). Having in 
mind Eq. (1) and the fact that the critical region (lines 1 and 2 in Fig. 1) along the 
X-axis is defined by |t(X)| <C 1, it becomes evident, that in the low-temperature range 
(X ~ X ), where [£ (X C )/A(T) < 1], the quantum subdomain gradually enlarges when 
the temperature is decreased or, equivalently, when X C (T) approaches X . In the zero- 
temperature limit (T — > 0) or, equivalently, for X C (T) — > X , this subdomain fills up 
totally the transition region. 

In the range of temperatures, where [£ (X C )/A(T)] > 1, the critical behaviour is totally 
classical. We should have in mind that the latter condition for a total classical criticality 
can be easily satisfied in the temperature range of quantum degeneration in real systems 
with a large zero-temperature correlation length (£ a). 

At high temperatures, where X c ~ 0, the scaling of correlation length £ can be conve- 
niently investigated by a scaling law with respect to X or (X — X c ), instead of t(X). The 
problems in the treatment of X-transitions in the high-temperature limit are analogous 
to those for the low-temperature T-transitions and the analysis can be carried out by 
following the ideas presented in Sec. 2.3. 

A quite special situation should exist when X = 0; see Fig. lb. In this case both T and 
X are equal to zero at the zero-temperature critical point, where the correlation length 
£(T, X) may exhibit scaling dependence on both t(T) and t(X). Obviously, such zero- 
temperature critical points will offer less opportunity for the quantum critical phenomena 
observation or, in some cases, this opportunity may not exist at all. To emphasize that 
the lack of any quantum critical phenomenon may be expected at such zero-temperature 
critical points, we have depicted in Fig. lb a classical transition region which completely 
fills up the total transition region (T, X) t , confined between the lines 1 and 2. 

However, there exist systems, where the critical temperature T c depends on X and tends 
to zero for X — > X = but the correlation length £ does not show any divergence with 
respect to the parameter X. In this case the function £(T, X), for all possible T and X, 
has the general form £[t(T), T C (X))], which describes only the scaling law with respect to 
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t(T) as given by Eq. (2) or Eq. (10). The transition region width along the T-axis tends 
to zero for T c — > and a transition region with respect to the parameter X does not exist 
at all. Thus the lines 1 and 2 in Fig. lb will definitely terminate at T = X = 0; see 
the dashed line 0-2. Certainly, the X-transitions at low temperatures offer much more 
favourable conditions for quantum critical phenomema, although the practical observation 
of such transitions in certain systems, where the parameter X cannot be gradually varied, 
is almost impossible. 

We end the discussion of the total and quantum transition regions by two remarks. In 
some systems at extremely low and zero temperature, the quantum transition region 
might be larger than the total classical region. Then quantum critical phenomena may 
happen outside the classical transition region. In this case, quantum critical phenomena 
may occur even if the size of the classical transition region tends to zero (the case depicted 
in Fig. lb). The second remark is about the possibility for a further generalization of this 
consideration by assuming that the symbol X represents more than one thermodynamic 
parameters: X = (X 1; ...). This generalization is straightforward. 

2. 5. Crossovers 

We have shown that the low-temperature critical behaviour can be either quantum or clas- 
sical depending on the specific (nonuniversal) properties of a particular system. Besides, 
there are no general arguments indicating that high-temperature and low temperature 
critical properties should be equivalent. Thus we may suppose that, in general, classical 
high-temperature classical low-temperature and quantum critical properties present three 
different types of critical behaviour. These three types of critical phenomena can be in- 
cluded in the framework of a quite general notion for a "high-low temperature crossover" 
(HLTC) of the critical behaviour which corresponds to the change of the critical properties 
(or of some of them) when the thermal length A varies from A < a at high temperatures 
to A > a (up to A > a) at sufficiently low temperatures, and vice versa. A special 
case of HLTC is the classical (g < 1) to quantum (g > 1) crossover (CQC) which de- 
scribes the difference between classical high-temperature phenomena and the quantum 
critical phenomena at T — > fTJ |20]- This CQC is treated by statistical methods 
(Sees. 4-7). 

The experiment has a finite accuracy and omits very narrow temperature intervals such 
as the asymptotic classical regions in the close vicinity of "almost-zero-temperature" crit- 
ical points (T c ~ 0). This should be taken into account in interpretations of experiments. 
Because of the zero temperature unattainability, quantum critical phenomena produced 
by zero-temperature X-transitions cannot be observed in a real experiment but certain 
low-temperature X— transitions might be experimentally investigated. On account of 
limitations in the accuracy of the equilibrium temperature measurement, the very nar- 
row classical region will remain unattainable for experimental studies and, hence, the 
experiment will reproduce results which, to some extent, will give an information about 
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asymptotic quantum critical phenomena at zero temperature. 
3. Ideal Bose gas 

3.1. Preliminary notes 

The model of IBG describes the phenomenon of Bose-Einstein condensation (BEC) jTJ 
ITT] ITU] ITT] . BEC of IBG is important for the undestanding of quantum phenomena in 
many areas of physics, in particular, in superfluid helium liquids excitonic phases in 
semiconductors [31] . the recently discovered BEC of dilute Bose gases in magneto-optical 
traps jH2 HHH HHl HUE EH] although in these and other real systems the interparticle 
interactions do affect the thermodynamic and correlation properties. 

We find reasonable to emphasize that the IBG does not contain interparticle interactions 
but it possesses two other properties which may cause a phase transition. Firstly, the 
quantum-statistical correlations act as attractive "pseudo-interactions." This can be seen 
from the negative sign of the first quantum correction to the equation of state of the ideal 
classical (Boltzman) gas coming from the Bose statistics [TUllllj . Secondly, the constraints 
imposed on the system usually have an effect similar to that of some interaction. Let us 
remember that the thermodynamics of IBG is ruled by constraints (a constant density, 
or, the less common constraint of a constant pressure and, why not, a constraint of 
constant temperature, which has been never investigated). These constraints are a simple 
experimental requirement and are used to define the chemical potential as a function of 
the temperature within the framework of the grand canonical ensemble [TU] ITT] . 

Besides, we have to emphasize that the ground states of IBG and NBG are quite different 
in their physical properties. While IBG produces the originally predicted by A. Einstein 
BEC, a very small interaction is needed to induce a superfluid ground state in NBG as 
demonstrated theoretically in a rigorous way by H. H. Bogoliubov j^Zj (see, also, Ref. 0). 

Another problem is the description of the phase transitions in Bose gases: from one side, 
the phase transition to BEC in IBG at finite and zero temperatures and, from the other 
side, the so-called A— transition, i.e., the phase transition to a superfluid state in NBG. 
These phase transitions exhibit quite different properties. Here we shall review the results 
from Refs. [23 ELU El S3 , where BEC in IBG has been investigated within the framework 
of the scaling theory of phase transitions [TJ H2] H3] . 

3.2. Thermodynamic equations 

The thermodynamic properties of rf-dimensional IBG, including BEC, are described by 
equations for the grand canonical potential Q and the number density p = (N/V) of 
bosons HUH]]: 




(12) 



k 
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and 

P = yJ2^ a k) ' ( 13 ) 

k 

respectively.Hete the brackets () denote a statistical averaging, (3 = 1/kgT, at and a^ 
are the creation and annihilation Bose operators for a plane-wave state of wave vector 
k = (hi, i = 1, . . . , d), V = (Li . . . L d ) is the volume of the gas. The self-energy e{k) = 
£o{k) +r is represented by the energy spectrum eo(k) = h 2 k 2 /2m of free (noninteracting) 
real particles or excitations, and the chemical potential \x = —r < 0; k = \k\. 

Throughout the paper we shall use periodic boundary conditions. The wave vector com- 
ponents ki = 2irli/Li, (li = 0, ±1,...), are given by the spatial dimensions Lj. They 
are supposed to be much larger than any characteristic length of the system, for exam- 
ple, Li ^> (£, A). In this case it is usually said that the dimensions Lj are "infinite". 
The stated condition allows to pass to the continuum limit, i.e., from a summation in 
Eqs. (12) and (13) to the corresponding integration. The quite uncommon case of com- 
posite Bose exitations due to long-range interactions in electron and magnetic systems 
can be included in the consideration by the generalization Eo(k) = h 2 k a /2m, (0 < o < 2) 
of energy spectrum e (k); see, e.g., Refs. [TJI13JI13]. ^ simple dimensional analysis of 
the exponent (3e(k) in Eq. (12) shows that £ = [h 2 /2mr) 1 ^ . 

With the help of the Bose distribution 

n(k) = (atar) = „ ,,\ > (14) 

Eq. (13) will take the form 

k 

The correct investigation of IBG thermodynamics in the continuum limit implies to take 
the density po — {N$/V) = n(0) of bosons with zero wave numbers out the sum (15). This 
is important for obtaining a correct description of the properties below finite temperature 
critical points T c > 0. For all other cases this separation is redundant. As we shall be 
mainly interested in the latter case, we shall avoid the mentioned separation; see Ref. [\2\ . 

Substituting the ^-summation in Eqs. (12) and (13) by a <i-dimensional integration, 

oo 

viWIW*"" (16) 

k 

K d = 2 1 - d /n d / 2 T(d/2), we obtain 

P = k B TA\- d g Wa+1) f-^j (17) 
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and 

-d. 



p = AX g (d/a) 



k B T 

Here P = — (fl/V) is the pressure, A = (2irh 2 /mksT) 1 ^ is the thermal wavelength with 
an exponent 9 = 1/cr, c.f. Eq. (1). The function g v (tz) and the parameter A which enter 
in above equations are given by the expressions, 



oo 



/ \ 1 /* oc doc , . 

*<(*) = W) j I > ( 19 ) 



and 

t77r d(i/2-i/^)r(d/2) ' 1 J 

(for a = 2, A = 1). 

The BE condensate is a coherent state of a macroscopic number No ~ N of bosons 
with a momentum hk = 0. The critical temperature T c of the transition to BEC in 
the momentum (k—) space is defined by the equation [r(T c )/T c ] = 0. The BEC order 
parameter O — ( a o ) / W is related to the square root of the number density p = (No/V) 
of the condensate bosons: \(f> \ = -Jpo- In the present problem, the ratio (3) takes the 
form 

Anr ^ 1/CT 



The IBG critical properties can be investigated with the help of Eqs. (17) and (18). They 
depend on the thermodynamic conditions imposed on IBG, i.e., on the way, in which the 
chemical potential fi = — r is determined. 

We shall briefly consider three constant pressure P jH] , a constant density p and 

spatial dimensions d > a [SHI HO] , and a constant density at dimensions d < a [12] ■ In all 
cases the critical regime is defined by the condition r < k B T and all critical phenomena 
are almost (r ~ k B T) or completely (r <C k B T) classical; c.f. Eq. (21). 

3. 3. Constant pressure 

The parameter r is calculated from Eq. (17) as a function of T and P. The result p(T, P) 
is substituted in Eq. (18) and, hence, one obtains the equation of state f(T,p,P) = 0. 
The solution T C (P) of Eq. (17) with r = yields the critical temperature 



T C (P) 



X d P 



C(d/a + l)Ak B _ 



<y/{d+a) 



[22] 



A characteristic feature of the phase transition at constant P is that the critical temper- 
ature is finite (T c > 0) for all dimensions d > and P > 0. The zero temperature critical 
point exists only in the limit P — > 0. This phase transition does not exhibit a crossover 
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to a low-temperature behaviour when the pressure P tends to zero because the critical 
exponents remain unchanged. There is a crossover, however, from Gaussian exponents at 
dimensions d > a to low-dimensional (d < a) Gaussian exponents. This crossover cannot 
be considered as HLTC or CQC. 

3.4- Constant density 

At a constant density p = (N/V), the equilibrium chemical potential n(T,p) is obtained 
as a solution of Eq. (18). The equation of state is given by Eq. (17) after the substitution 
of solution r(T,p). Eq. (18) with r = yields 



Tc(p) 



(23) 



and, therefore, T c (p) > for d > a, provided p > 0, whereas T c = for d < a; note, that 
the zeta function ((d/a) tends to infinity when the ratio d/a decreases to unity. 

The finite temperature critical behaviour (T c > 0) is identical to that of the classical 
Berlin-Kac j3H] spherical model [SHIED! (for the spherical model see, e.g., Ref. [33], and 
the brief note in Sec. 4.1). The reason for this behaviour is in the constant density p 
condition which, as is seen from Eq. (13), is equivalent to the(mean) spherical constraint 
on the variations of mean densities {a^a^} /V. The critical exponents are shown in Table 
1, including the case d > la when they take Gaussian values [HE!. All critical exponents 
are given by their usual notations and definitions; see, e.g., Refs. [T2| IT5]. 

Table 1. Values of critical exponents (p = const). 
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The critical exponent 7 describes the off-diagonal susceptibility 

* = -(^L~ (i/r) - (24) 

where h is a complex number - a fictitious external field conjugate to the Bose operator 
(^2% at). For T c = 0, a new critical exponent 7 can be introduced to describe the density 
n(0) ~ ksTx- For d < a, the values of 7 from Table 1 and the relation 7 = (7 — 1) yield: 
7 = d(o~ — d) for d < o~, and 7 = 00 for d = o. The exponent 7 does not exist for T c > 0, 
where the density n(0) ~ x is described by the usual exponent 7. 

In Table 1, the exponent a s is an auxiliary exponent intended to give a more detailed 
description of IBG specific heat for d > a. For d > a, the specific heat C(T) of IBG and 
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spherical model can be represented in the form 

C(T) = C + C(T - T L )~ as , (25) 

where Co is a regular constant. This behaviour corresponds to the zero value of the 
usual specific heat exponent a, defined by the scaling law C(T) = C^\T — T c \~ a , where 
C(°) is the scaling amplitude. Usually several shapes of C(T) are usually ascribed to the 
value a = 0, namely, a logarithmic divergence (ln|T c — T\/T c ), a finite jump as is in the 
MF theory, and cusps of different forms. In order to avoid this arbitrariness, Fisher [13] 
introduced the second term in Eq. (25). The values of a s in the Table 1 are negative 
and this next-to-leading term in C(T) is always small but the derivative dC(T)/dT is 
divergent at T c and this is an information about the cusp shape at T c . 

The essential difference between the properties of the phase transition to BEC in IBG 
and the usual second order phase transitions can be easily seen by the comparison of the 
Landau free energy PQIH] F(<p) = Vf(ip) of a standard second order phase transition, / = 
(ro(p 2 +ip 4 ) with the free energy density f{ip) of the Bose gas at constant density jlQ], given 
by / = (r + ip 2 ) 3 , where r ~ t(T). The simple form of these free energies corresponds 
to the suitable choice of units for / and the order parameter ip. The functions f(<f) are 
shown in Figs. 2a and 2b for several values of the parameter r . While the well defined 
minima of / in Fig. 2a describe the ordered phase of usual second order phase transitions 
BEC is described by the inflection points (df/dtp = d 2 f/d^ 2 = 0, d 3 f/dip 3 > 0) of 
the /—curves shown in Fig. (2b). At these inflection points the /—curves intersect the 
ip— axis; the regions where / < 0, see Fig. (2b), correspond to a positive value of the 
chemical potential and, hence, are unphysical |30] . 

Certainly, the properties of IBG below the condensation point r < also correspond 
to a continuous phase transition which is quite different from the standard second order 
transition. Note, that the coefficient of the <p 4 — term in the free energy density of the 
Bose gas, / = (r + <^ 2 ) 3 ? is proportional to t and tends to zero when t — > 0. This 
behaviour of the fourth-order term is quite similar to that of the Landau expansion for 
tricritical phenomena [1J. But the behaviour of IBG differs from that at standard tricritical 
points PQ, where the temperature dependence of the Landau coefficients is different. 

3.5. Zero temperature condensation 

The zero-temperature (T c = 0) BEC at a constant density was investigated in Ref. jl2*] . 
This condensation is possible at a finite constant density (p > 0) and a low spatial 
dimensionality d < a. The results for the critical exponents at d — a and < d < o are 
shown in Table 1. These critical exponents describe the scaling laws with respect to the 
variations of temperature, for example, C(T) = C /T a . 

It can be seen from Table 1 that the zero-temperature condensation at a constant density 
exhibits quite unusual critical properties. The critical exponents 7 and v for d < o can 
be obtained from the familiar values at a < d < 2a by the change of sign of (d — a). The 
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Figure 2: (a) The function f((p) for a standard second order phase transition for: ro = 1.0 
(curve 1), r = 0.5 (2), r = —0.5 (3), r = —1.0 (4). (b) The function f(<p) for the 
IBG at constant density (the curves 1-4 correspond to the same values of r as given for 
Fig.2a). 

Fisher exponent rj has its usual value but the exponent a is negative rather than zero as 
is for d > a. For zero-temperature T-driven transitions (Sec. 2.3), as is in our case of 
d < cr spatial dimensions, the negative values of the critical exponent a are consistent 
with the Nernst theorem. When d > a, a = and this theorem is satisfied in the limit 
T c (p) — > of an extreme dilution (p ~ 0) because the scaling amplitude Cq tends to zero. 

The infinite values of the exponents v and 7 at d = a indicate the exponential divergence 
of the correlation length £ = (ft 2 /2mr) 1//<J and the susceptibility x = (1/r) at T c = 0. 
These exponential divergences are known from the mechanism of critical fluctuations in 
classical systems at their lower borderline (critical) dimensionality d^. Here the exponen- 
tial divergence of £ and \ comes from the same low dimensional effect. The exponents 
given in the last two lines of Table 1 can be referred to as low dimensional spherical 
exponents. 

The next question is why this zero-temperature critical behaviour exists. It is intuitively 
clear that the low dimensional (d < a) IBG should have a ground state of a total fc-space 
condensation (iVo = N) at T = which means that the solution T c = of the equation 
[r(T c )/T c ] = should exist. For T ~ 0, (r/ksT) 1 and, hence, the contributions to the 
sum (15) are given by terms with low energies e(k) ~ r, i.e., with small wave numbers k. 

There is an onset of BEC at T ~ 0, where the bigger part of bosons are at states 
e(k) ~ but the macroscopic condensation occurs only at T = [42J. At T = the order 
parameter 0o jumps from zero to |0o I = p 1/2 and all N bosons enter in the ground state 
(k = 0). This peculiar critical behaviour results from the constant density p condition 
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(the spherical constraint). 

The same behaviour can be obtained from the classical spherical model. Although the 
fc-space condensation is a pure quantum effect due to the Bose statistics, the critical 
behaviour around critical points (T c > 0) is totally ruled by classical fluctuation effects 
described by the classical spherical model. 

The IBG critical behaviour can be exactly derived from thermodynamic Eqs. (12) and (15) 
because this model is exactly solvable. Sometimes the RG analysis of such exactly solvable 
models is also useful, for example, in studies of the crossover of the critical behaviour of 
IBG to that of NBG [23] . 

Note that a zero-temperature BEC (T c = 0) is possible also in the limit p —>■ (extreme 
dilution) for the case of constant density, and in the limit P — > for the case of constant 
pressure for each spatial dimensionality D > 0. As mentioned in Sec. 3.3 at these onset 
transitions there is no HLTC or CQC. 

3.6. Classical and quantum critical regimes 

The critical behaviour of IBG in a close vicinity (r <C k B T) of phase transition points 
{T c > 0) to BEC of IBG is definitely classical. All our calculations have been performed 
by keeping the leading power in {r/ksT) in the corresponding series for Bose integrals (19) 
in powers of {rjk B T) < 1. This approximation has been used to determine the leading 
scaling behaviour, namely, the scaling laws. This essentially classical expansion gives the 
critical regime {p ~ 0). 

The transition region (Sec. 2.2) where the phase transition phenomena occur, can be 
defined for a broad interval of temperatures around T c > 0, i.e., by r < k B T '. In this 
region but not very near (r <C k B T) to T c we must take into account some next-to- 
leading terms in the same series expansion of the Bose integral (19) in powers of r/k B T. 
These secondary terms have a quantum origin and, therefore, taking several of them as 
corrections to the leading scaling powers, we shall obtain the quantum corrections to the 
main scaling behaviour. 

Thus the classical region (a-0-b) discussed in Sec. 2 can be defined by the inequality 
r <C k B T. For both T C {P) and T c {p) this region will look like the shaded wing 1 — — 2 
in Fig. lb; see the dashed line — 2. The quantum region will be outside this vicinity of 
critical line, in the rest part of the critical domain (r < k B T) discussed in Sec. 2. 

The parameters P and p rule the critical temperature but they do not participate in 
the scaling laws. Thus they have not a direct effect on the critical region size. This is 
a particular IBG property. The size of the classical critical region is reduced with the 
decrease of P or p, i.e., of T c . This is valid for both widths of classical (r <C k B T) and 
total (r < k B T) critical regions along the X— axis; X = {P,p). In fact, at fixed X we 
should consider the inequalities r < k B T c {X) and r <C k B T c {X). When T C {X) is lowered 
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the critical regions width along the X-axis will decrease to zero for X — > 0. 

The variations of the overall density p in the ground state [T = T c (p) = 0] of the low- 
dimensional (d < a) IBG does not produce critical effects, because the condensate density 
is equal to p (all particles in IBG are in the condensate). If the zero-temperature critical 
point, containing the complete BE condensate of IBG is approached by T-variations in 
the disordered phase, the observed critical phenomena will be mainly classical on account 
of the real phase transition to BEC at T c = 0. 

The phase diagram of the zero-temperature BEC in low dimensional systems at a constant 
density (Sec. 3.4) is very simple because the critical line coincides with the p-axis. The 
critical line [T c (p) = 0, < p < oo] can be approached only by T-variations and the 
critical phenomena are purely classical. 

Finally we shall note that the criterion (9) is not satisfied for the above mentioned tran- 
sitions. There is no real CQC in IBG because of the lack of fluctuation interactions but 
we can consider the crossover from the high-temperature to the zero-temperature at a 
constant density p > caused by the change of the dimensionality from d > a to d < a. 
This formal dependence describes the difference between the high-dimensional (d > a) 
critical phenomena at finite critical temperatures and those of zero-temperature BEC at 
a constant density p > and low dimensionality (d < a). 

4. Interacting bosons 

4-1. General methods 

The interacting bosons can be investigated by propagator expansions within the Green 
functions method [HI ESJ EZ1 EH] • The infinite perturbation series are truncated on the basis 
of several approximations. The perturbative approach breaks down in a close (asymp- 
totic) vicinity of the critical point because of the strong fluctuation interactions. In 
this situation, after accomplishing the MF analysis, the strong fluctuation interactions in 
the Ginzburg fluctuation region PJ around the critical point and their effects on the 
asymptotic critical behaviour are investigated by RG [T| \V2\ IT3]. 

The RG method consists of a suitable choice of a length-scale transformation combined 
with the so-called loop expansion (see also Sec. 4.2). In the framework of the general 

RG approach the MF approximation is equivalent to the so-called "tree approximation" - 
the lowest order theory within the loop expansion. Certain aims as, for example, the eval- 
uation of the size of the fluctuation effects by a calculation of the Ginzburg critical region 
size or investigations of qualitative features of fluctuation effects in the pre-asymptotic 
transition region, can be achieved by suitable perturbation calculations without RG (see, 
e.g., Ref. [13 HE], and Sec. 6.3). This outline of general investigation methods of critical 
phenomena is common for both classical and quantum statistical models. 

For the interacting bosons we must use the Hamiltonian H[ip(r)} in terms of second- 
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quantized Bose field operators ip(r) which for practical calculations are expanded in terms 
of annihilation and creation operators mentioned in Sec. 3.1. Within this formulation 
which gives the opportunity to use the powerful Green functions method the quantum 
effects are "hidden" in the commutation relations for the field operators ip{r) and, hence, in 
the time ordering procedure under the trace for partition function and statistical averages. 
For finite temperatures the "time" ordering procedure, i.e., the account of quantum effects, 
is performed with the help of an auxiliary "time" r, reciprocal to the temperature: (0 < 
r < hf3). This quantity is called an "imaginary time" because of the relation r = it with 
the real time t. 

Alternatively, the theory can be formulated by the Feynman path integrals I5TH loT} 152*] . 
With the help of the coherent state representation the grand canonical partition function 
can be written as a functional integral over a (c— number) field <f)(x) depending on a 
(d + 1) -dimensional vector x = (t, r) in restricted "time"-space: [(0 < r < h/3),r e 
V]. So, we may introduce a (d + l)-dimensional "volume" V^+i = (/3Li...L d ), where 
f3{= L Q ) is the finite size ("thickness") of the time-space "hyperslab". Therefore, the 
quantum effects bring an extra dimension hf3 along the new (r— ) axis. At sufficiently 
high temperature (classical limit), the r— dimension collapses to zero and the behaviour 
of the system is classical. Of course, the same classical limit can be achieved by setting 
h = 0. At sufficiently low temperatures the r— dimension is large and becomes infinite 
at the absolute zero (T = 0). At such temperatures the classical limit (h — > 0) does not 
produce reliable results. 

The field <p(x) represents the order parameter, which describes the ordering below the 
phase transition temperature. The equilibrium ordering is given by the statistical average 
< <p(x) >. The fluctuation phenomena are represented by the fluctuation part S(j)(x) = 
[4>{x)— < 0(x) >] of the (nonequilibrium) field 0(x). In usual cases, for example, when the 
phase transition is a result of a spontaneous breaking of the global discrete or continuous 
symmetry, the average < <f)(x) > does not depend on x, i.e. we have a uniform equilibrium 
order parameter. For the sake of simplicity we shall often use this example. 

The complex Bose field <p(x) is sometimes referred to as a "classical field". This term 
points to the difference from the field operators which obey commutation rules. In fact, 
the classical fluctuation field is 0(0, r) = 0(f) and corresponds to the classical limit 
(h = 0). The r-variations of the field <fr(x) are created by quantum fluctutions. The way, 
in which the quantum fluctuations take part in this picture is very similar to that, in 
which the classical fluctuations appear as spatial variations of the field 0(x). 

Furthermore, in order to describe a large variety of systems one may consider 0(x) as 
a (n/2)— component complex vector field: <f)(x) = {<f) a (x),a = 1, ...,n/2}, where the 
compoments <p a are complex functions of x and n is an even positive number depending 
on the symmetry of the system, or, more precisely, on the symmetry of the possible ground 
states (orderings). The number n is often called a "symmetry index" or a "number of 
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order parameter components." 

For certain problems the field theories of n— component real fields are equivalent 

to the theories of (n/2)-component complex fields. The classical systems are usually 
described by real scalar or vector fields whereas the quantum systems are represented 
either by real field components or by complex field components of twice less number 
(the choice depends on the aims of the particular investigation). Thus one may compare 
the results for quantum systems with a symmetry index of ordering n/2 to the results 
for classical systems with a symmetry index n. Moreover, this correspondence can be 
extended to all integer numbers n > 0. Then a complex field with n' = n/2 components 
will correspond to a classical Ising system for n — 1, i.e., n' = 1/2. A Heisenberg 
ferromagnet (n = 3) will be described by a complex vector field of n' = 3/2 components. 

On account of the analytical relation r = it the purely quantum fluctuations reveal the 
intrinsic quantum dynamics of the system. The dependence of physical quantities on real 
time t is obtained by an analytical continuation (r — > it) of the results from calculations 
to the real time axis t. While in classical systems the static and dynamic phase tran- 
sition properties are considered separately, in the framework of the quantum statistical 
mechanics these phenomena have a unified treatment. The intrinsic quantum dynamics 
contained in the quantum statistical correlations is important for the dynamical critical 
properties of Bose systems; see, e.g., Refs. j53J EU E3 EH| • Note, that this intrinsic quan- 
tum dynamics can be substantially changed under the action of external time-dependent 
potential(s) [HJI22I- 

Adopting periodic boundary conditions along the r-axis, 0(r, r) = <p(r + hj3,r), and 
having in mind the periodic boundary conditions along the spatial axes (Sec. 3.2) we can 
expand the field components 4> a (x) in Fourier series 



classical fluctuations only. In the x-representation this corresponds to the limit (3 — > 0, 
in which the system size ("thickness") along the r-direction tends to zero. 

The functional formulation of Bose systems in terms of (c— number) complex functions 
<p(x) is performed directly for the original microscopic second-quantized Hamiltonian be- 
cause the commutation relations of the Bose operators ip(r) have a direct classical limit 
O HU EH] . The microscopic Hamiltonians of fermionic and spin systems are transformed 
to effective Bose Hamiltonians with the help of Hubbard-Stratonovich transformations 
and Feynman path integration EUl EDI OHJ E2J ■ The effective field Bose Hamiltonians 
have the same meaning and role in the theory of quantum phase transitions as the effective 




(26) 



where the (d + l)-dimensional vector q = (ui, k) is given by the Matsubara frequencies 
uj\ = 2nlkBT/h, where / = (0, ±1,...); qx = (uit + k.r). Setting u\ — in Eq. (26) we 
neglect the quantum fluctuations and the amplitudes [<fi(0,k)/\f/3 — > 4>(k)) will describe 
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field Hamiltonians for classical systems EE] • In both cases of classical and quantum 

systems the effective field Hamiltonians is not exact counterpart of the respective original 
Hamiltonians but they provide a correct description of the (quasi) macroscopic phenom- 
ena which is enough for a reliable treatment of phase transition problems. The lack of 
an entire correspondence with the microscopic model in nontrivial cases of interacting 
particles is a result of certain forms of coarse graining which are unavoidable product of 
the known theoretical techniques of derivation of effective field models. 

As our interest is focussed on standard second order phase transitions we shall consider 
(/> 4 -effective Bose "Hamiltonians" TC of fermionic and spin systems, as well as, the same 
type microscopic Hamiltonians of genuine Bose systems. In fact, the "Hamiltonian" 7i 
that we shall investigate is related to the action S of the system by 7i = — <S; here we 
follow a conventional terminology used in a number of papers; see, e.g., Ref. [20| . 

A quite general quantum effective Hamiltonian {H = (3H) comprising a number of systems 
can be written in the form 

= E^o 1 ^)!^)! 2 + ^ E C(?i)^(92)0a(g 3 )^(gi + ? 2 -g3), (27) 

where v > is the interaction constant. The (bare) correlation (Green) function Go(q) = 
(I0a(<?)| 2 )o is given by 

G \q) = iun+eik) (28) 

for interacting real bosons and transverse XY model [HJ For other effective Bose 
models of quantum systems the bare correlation function Go(q) can be written in the 
form HH 

I. , Im 

G o\Q) = ~^r + cka + r (29) 

with positive exponents m,m', < o < 2, and Landau parameters c > and r ~ 
t(T) or r ~ t(X). The parameter c can be always presented as c = h 2 /2m(g) where 
the effective mass m(g) of the composite bosons (Bose excitations) depends on some 
interaction constant g. 

The microscopically formulated Bose Hamiltonians contain an upper cutoff A = (ir/a) 
corresponding to the Brillouin zone ends (—n/a) < ki < (n/a). An example of such 
Hamiltonians is considered in Sec. 5. The effective quasimacroscopic Hamiltonians derived 
from microscopic Hamiltonians contain only large-scale spatial fluctuations of the field 
4>(q) and, hence, the corresponding cutoff A is relatively small: < A (tt/&) [45J. 

An example of a quasimacroscopic (effective field) Hamiltonian, derived by a Hubbard- 
Stratonivich transformation from the microscopic TIM Hamiltonian, is discussed in Sec. 6. 
As far as we are discussing critical phenomena (£ 3> a, A ^> a) the precise value of the 
cutoff A is not important. It is, however, important to have the conditions £ > (1/A) and 
A > (1/A) satisfied and then, the relevant phenomena included into consideration. 
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The statistical treatment implies the calculation of the generalized grand canonical parti- 
tion function 

Z(T, a) = Jv<j) e - H ^ ; (30) 
which gives the Gibbs free energy 

Q = -p\nZ(T, a) . (31) 

In Eq. (30), a denotes the Hamiltonian parameters (c,r,v), J D<p denotes the functional 
integration over all allowed field (j)(x) configurations 

III II dfc(x)d4> a (x) , (32) 

a=l zeV (d+1) 

or, in the fc-space, 

/ Ud^MdMl) (33) 

a;q 

over all allowed Fourier amplitudes 4> a (q). Constraints on the field configurations lead to 
a critical behaviour change. Within the present functional formulation the constraint of 
a constant density (13) is given by the mean spherical condition 

^ = £<i^(?)r> • (34) 



a.q 



Therefore, the mean square of the Nn/2 dimensional complex vector <p = [4> a (q)] (with 
components given by all possible a and q) describes a sphere of radius Nn/2. In the 4 - 
Hamiltonian of the nonideal Bose gas (NBG) given by Eqs. (27) and (28) with r = — /i, 
this constraint can be taken into account in the so-called large-n limit (n oo) [T| IT2*]. 

Let us note, that for certain problems JoTj the effect of constraints like (34) can be taken 
into account in Eq. (27) by adding an auxiliary 4 -interaction term: ~ m0 4 . In contrast 
to the ^-interaction terms 4>a( ( li)4 > ^(q2)4>a(q3)4 ) i3(qi + <?2 — 93) in Eq. (34), the auxiliary 
interaction is given by two summation g-vectors: 4'* a {qi)(t>*p{q2)4'a{qi) ( t ) p{.q2) |57J. 

The practical calculations are carried out by a substitution of the summation over the wave 
vector components k{ with an integration as shown by Eq. (16), provided the correspond- 
ing dimensions Li satisfy the criterion Lj 3> £ for a quasi-infiniteness. The summation 
over the Matsubara frequencies oji can be substituted with the integration 



00 

1 f duj . 

iZ-lwv (35) 

1 —00 

provided the temperature T is low enough or, in an exact mathematical sense, if T — >• 0. 
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In order to determine the condition, under which the integration (35) can be used without 
substantial errors in the final results, let us consider the model (27) corresponding to 
the correlation function Go(q) given by Eq. (29) with w! — and m — 1. Then the 
^-dependent modes <p(ui ^ 0, k) in Eq. (27) will yield relevant contributions to the 
partition function, if they have a relatively high statistical weight. This may happen 
provided U\ = lnk^T < r near the critical point r ~ 0. In terms of characteristic lengths 
A = (47rc//cB^) 1//<J , and £ = (c/r) 1 /' 7 we have A > {Sti 2 ) 1 ^^ which is consistent with the 
criterion (3); here we have used Eq. (2) and c = h 2 /2m. 

If the criterion (3) is satisfied, we can substitute the summation over ui with an integration 
as shown by the rule (35) without introducing a substantial error in the calculations. If 
the criterion (3) is not satisfied, the frequencies u>i can be neglected and the quantum 
fluctuations ignored. 

There exists a well developed theory [SDl EE EH HUB EH E2] of the A— transition in 4 He 
which can be applied to other Bose fluids with critical temperatures T c > 0. We shall 
not dwell on this important topic. Our attention will be concentrated on the case T c — > 
when the asymptotic scaling is influenced by quantum effects. 

4-2. Notes about renormalization group 

A number of problems in the scope of this article are usually solved by RG in the q— space, 
i.e., by field theoretical variants of RG [TJ H21 EH EE] . However, because of the simulta- 
neous presence of more than one effect (quantum and thermal fluctuations, disorder, 
anisotropy or gauge-field effects) a comprehensive RG investigation of quantum systems 
can be performed by widely applicable and not extremely sophisticated variants of RG 
such as, for example, the e— or the (1/n)— expansion PQ^JEI]- The latter expansion 
is convenient for the calculation of critical exponents but is quite hard as a method of 
investigation of the RG transformation [T2*] . 

In the e = (du — d) —expansion the difference (du — d) is used as a small paremeter; d\j is 
the so-called upper critical (borderline) dimensionality, above which the fluctuations are 
irrelevant and the system is described by the MF approximation. For the classical variant 
(uji = 0) of the Hamiltonian (27) we have du = 2a but for other Hamiltonians du may 
have another value pQ. Within the e— expansion the physical quantities can be calculated 
as series in powers of e and the order of accuracy of the calculation (first, second, etc. 
order in e) coincides with the order of the loop expansion which has been taken into 
account [H CHj , that is, the calculated order in e exactly corresponds to the respective 
order in the loop expansion. The method works for dimensionalities d > d^, where du 
is the so-called lower critical (borderline) dimensionality, below which the ordering is 
destroyed by the fluctuation effects; for the Hamiltonian (27) with (a;; = 0) we have 
d L = a p. 

The e— expansion has been widely used in investigations of critical phenomena and in the 



26 



remainder of this paper we shall discuss results obtained by this method. Note, that the 
e— expansion is asymptotic and all reliable predictions about the possible types of fixed 
points (FPs) of the RG differential equations [TJ H3] , or, equivalently, of the Wilson-Fisher 
recursion relations [TJ H2J EH1 EH EH3 EH] are obtained within the one-loop approximation 
(equivalent to the first order in the e— expansion). 

The next order of the theory, i.e., the two-loop order, which is equivalent to the second 
order in the e— expansion, provides better quantitative results for the critical exponents, 
corresponding to the FPs, but this and even the higher orders in the loop expansion do 
not essentially contribute to the understanding of the main qualitative features of the 
system, such as the possible types of critical behaviour and their stability properties. The 
same is valid for the precise determination of the domains of attraction of simultaneously 
existing stable FPs of complex Hamiltonians describing more than one interaction effect; 
this point has been discussed in Ref. |HZj on the basis of a quite complex model. 

The above notes justify our approach to the classical and quantum critical phenomena in 
complex systems with competing effects, for which we shall often use results in first order 
of e. We would like to emphasize an important but not widely accepted feature of the 
modern RG methods, namely, that the lowest (one-loop) results reveal the entire picture 
of the possible critical phenomena whereas the higher-order considerations in the loop 
expansion are useful for other aims: the improvement of numerical results for the critical 
exponents and also the investigation of the asymptotic nature of the RG perturbation-like 
method(s). 

Bearing in mind that the most often used expansions in the theoretical physics are asymp- 
totic we should not take the asymptotic nature of the e— expansion as a great disadvantage, 
or, as a reason for an unreliability of the e— results. The experience accumulated by the 
research done during the last 20 years firmly indicates that the e— expansion results are 
meaningful and quite useful, provided they are consistently analyzed and interpreted |[7] . 

The main types of critical behaviour which can be described by a concrete model (Hamil- 
tonian) are given by the FPs of the respective RG equations. Different FPs describe 
different types of critical behaviour and, hence, different critical exponents. These expo- 
nents are obtained from the relevant and irrelevant stability exponents of the FPs. If a 
FP appears of order 0(e 2 ) but does not exist within the one-loop (first order in e) RG 
equations, it cannot be accepted as a new reliable object of investigation and the reason 
lies in the asymptotic nature of the e— expansion. 

At relatively high spatial dimensionalities d > du which usually do not correspond to real 
systems the RG equations have one stable FP, the so-called Gaussian FP (GFP). This FP 
describes the MF behaviour for d > dy. The same GFP is unstable for d < du and the 
stable critical behaviour in this domain of spatial dimensionalities, which includes the real 
system dimensionalities (d = 1, 2, 3), is described by another stable FP. This is usually a 
nontrivial FP which gives e— corrections to the MF values of the critical exponents, i.e., 
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this FP yields a nontrivial (non-MF) critical behaviour. Sometimes such FP is called 
the Wilson FP EU EE] but the most often used terms are either Heisenberg FP (for 
systems with a continuous symmetry, n > 1) or Ising FP (for systems with a discrete 
symmetry, n = 1). When the system has a XY symmetry (n = 2) the nontrivial FP is 
called XY FP. GFP has zero coordinates (ra = vq = 0) in the (r, v) parameter space 
of the Hamiltonian but the nontrivial FPs have nonzero coordinates which are given to 
first or higher order in e, depending on the order of the loop approximation used in the 
particular investigation. 

The possible types of critical behaviour predicted by RG are classified in universality 
classes. For a given structure of the Hamiltonian, one can define the universality classes 
by the couple (d,n), that is, the critical behaviour depends on the dimensionality d < djj 
and the symmetry index (number of components of the order parameter field <p) n. But 
one must keep in mind that the change of the mathematical structure of the Hamiltonian 
usually leads to another double series (d, n) of universality classes. The equivalence of 
the critical behaviour of many substances, described by the same Hamiltonian structure 
and the couple (d, n) is recognized as a property of the universality of critical phenomena. 
The MF behaviour at d > djj is superuniversal because it does not depend on d and n. 

Below the upper critical dimensionality djj more than one FPs may exist and be stable in 
different domains (domains of attraction) of the Hamiltonian parameter space. Usually 
this happens for complex systems with competing effects. The Hamiltonians of such 
systems include more than one interaction term. The simultaneous stability of more 
than one FP is a manifestation of the simultaneous action of more than one fluctuation 
interaction or the presence of other effects which interfere with the main fluctuation 
interaction; the latter is usually represented by a 4 — term in the Hamiltonian as shown 
in Eq. (27). This situation is quite common in the RG analysis of Hamiltonians with 
competing effects. Now one is faced with the hard task to determine precisely the FPs 
domains of attraction which is often necessary for the correct prediction of the critical 
behaviour. We have already mentioned that for quite important cases this cannot be 
reliably done even in orders higher than the first order in e jHZj. In other cases, the effects 
competing the main fluctuation interaction ((f) 4 ) are hidden in additional modes like the 
modes <p(uJi, k) with ui ^ of the quantum fluctuations, and under certain circumstances, 
these modes may have an essential effect on the critical behaviour. 

The RG predictions are uncertain when there are no stable FPs. This may happen for 
a definite domain of the Hamiltonian parameter space or, in extreme cases (see Sec. 7), 
for the whole parameter space. The lack of stable FP is usually interpreted as a lack of a 
stable (multi) critical behaviour of usual type. Sometimes, the absence of stable FP can 
be considered either as a signal for a first order phase transition or as a more particular 
type of a continuous phase transiton. The reliable prediction of the phase transition type 
in this situation requires additional heuristic arguments and investigations outside the 
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scope of the RG method; see, e.g., Refs. [T] IT21 ITBI ITo] . 

4-3. Classical-to- quantum dimensional crossover 

We have already mentioned that the <i-dimensional quantum systems resemble D = 
(d + l)-dimensional classical ones. This is valid for the model (29) with m! — and 
m = a. By the formal substitution k^T = c 1 ^ / L in cui, where L = A = (c 1 ^ /kgT) 
is the thermal wavelength corresponding to the model (29), we obtain a D = (d + 1)- 
dimensional momentum q = [qf,i = 0, l,...,d] with components qi = (2irli/Li). The 
representation gives a total formal correspondence between the present time-space geom- 
etry and a hyperslab with (d + 1) spatial dimensions. This formal analogy was revealed 
for the first time in Refs. [HHj. H2I] and further explored in Ref. |70j . Using the formal 
analogy we have shown that for this effective model the notation (1) yields the exponent 
9 = 1. 

The temperature dependence A = Ao/T 1//fT , as in IBG, is valid for other effective models, 
for example, for the model (28) and that given by Eq. (29) with ml = and m = 1. 
Therefore, the form of the scaling law A = Xo/T 6 depends on the values of exponents 
m, ml and a: 9 = 8(m,m' ,o~). A simple dimensional analysis of Eqs. (28) and (29) 
shows that 9 = 1/z. The same result (9 = 1/z) can be obtained by the RG rescaling 
transformation pJElEj- Then the criterion (9) becomes zuq < 1. Usually the (bare) 
values of exponents 9 and v = 1/er of A and £, respectively, receive perturbation corrections 
which are calculated with the help of RG. 

It has been shown by Hertz [20] with the help of the above mentioned method as well by 
a RG analysis that 

z = ^ (36) 
m 

and, hence, that the real CQC is given by 

D = d + z. (37) 

CQC is always associated with the dimensional crossover (37) which was revealed before 
Hertz by Pfeuty and Elliott |T7j and by Young ^J] on the basis of TIM (for the latter, 
m' = 0, a = 2, m = 2, and 2=1). Note, that the upper borderline dimensionality 
changes from its classical values du to the "quantum" value d^ = (djj — z) and the lower 
borderline dimensionality becomes d° L = max[0, dt — z)]. 

Thus we should expect that the quantum critical phenomena in d-dimensional (quantum) 
systems are described by the universality class (D,n) known for D = (d + z)-dimensional 
classical systems. If the shift d — ► D of spatial dimensionality d is the only result of 
quantum effects, we can say that the quantum critical phenomena in the particular d- 
dimensional quantum system (or a class of systems) obey a form of universality, namely, 
that they are identical to the critical phenomena in the corresponding D-dimensional 
classical system (or class of systems). Under the term "corresponding classical system" 
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of the quantum model (27) we understand the system described by the classical variant 
(wj = 0) of Eq. (27). 

This definition of the universality of quantum critical phenomena allows a comparison 
between the critical properties of classical and quantum systems. There are other points 
of view on the definition of the universality for quantum critical phenomena. For example, 
sometimes it is said that the quantum critical phenomena in a system are universal, if 
a classical critical behaviour can be found at the same or another spatial dimensionality 
which describes the same critical phenomena [HJ E21 EH1 • This point is important, for 
example, in the discussion of the extraordinary zero-temperature critical behaviour of 
NBG (Sec. 5). 

5. Nonideal Bose gas 

5.1. Preliminary notes 

The first studies [JJJ [JJ3 [THJ E7J of NBG by RG were performed in the framework of the 
operator formalism El IHH E7J ■ The direct RG applications to the second-quantized 
NBG Hamiltonian lead to a scaling law for the commutation relations of the field operators 
V>(r). In Refs. [ZU UJ3 [23 EJJ this was interpreted as a rescaling of the boson mass m 
which grows with the successive RG transformations. The fact that the quantity, to 
which such a rescaling should be ascribed is the thermal length A or, equivalently, the 
factor (1/mT) associated with it, has become clear later. However, this circumstance 
was undoubtedly irrelevant to the main conclusion of these works, namely, that the finite 
temperature (T c > 0) BEC should exhibit a universal critical behaviour corresponding to 
the universality class (d, 2) of classical XY model. The latter result has been confirmed 
and extensively discussed in other studies; see e.g., [?H1 E3 EDI EH IH3 E3 EH E3 An 
error in the calculations [Til [73] which is irrelevant for finite critical points (T c > 0) but 
very important for the zero-temperature critical behaviour at T c = 0) has been pointed 
out in Ref. [HI] and noticed also in Ref. [23J. Other relatively early papers on critical 
properties of real bosons systems are mentioned in Refs. [Tj IH| l2~3~|. 

The zero-temperature critical behaviour of NBG was treated for the first time in Ref. [HE] 
in the one-loop approximation and despite of the correct result for the critical exponents 
there are errors in the derivation of the RG equations and their FPs. In Ref. jHZj the 
correct calculation of the self-energy function, the Fisher exponent 77 and the dynamical 
critical exponent z was made in the two-loop approximation. 

The Wilson-Fisher recursion relations PIE] appropriate for a correct and thorough treat- 
ment of the critical properties of NBG at low and zero temperature critical points have 
been derived in Ref. jHIj and in Sees. 5.2-5.6 we shall follow this approach. The per- 
turbation series and the diagrammatic representation of perturbation terms are standard 
and we shall not dwell on technical details; for a more instructive explanation of the 
corresponding calculations, see Refs. [H EH EH] • 
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Because of the model particular properties the low-temperature limit gives an extraordi- 
nary opportunity for an exact summation in all orders of the loop expansion and, there- 
fore, for obtaining of the low-temperature critical behaviour without approximations. 
This was demonstrated for the first time in Ref. |81j . However, the exact result [HI] has 
been overlooked in subsequent papers [23 EE EH1 OH] where the same low-temperature 
critical behaviour was rediscovered in the one-loop approximation and applied to dilute 
Bose systems. In Ref. [HE], the Gaussian-like critical behaviour [HI] to one-loop order 
(see also Sec. 5.3) was very conveniently interpreted as a form of "quasi-universality." In 
Ref. [23*] I9~T] the RG recursion relations [HI] were re-established in the form of one-loop 
RG differential equations and solved by the Rudnick-Nelson [92J method of intergation. 

More recently, the simple structure of the perturbation series discussed below, was found 
also in the nonrelativistic scalar field theory and appropriately called the "scaling anomaly" 
[HH]- Sachdev and co-workers [94 J introduced another appropriate term - "zero scale- 
factor universality" - for the same exact solution [HI]. These authors overlooked the 
opportunity to take advantage of perturbation series exact summability (see below) but 
recently the exact result has been rediscovered and lengthy explained by one of them 
without any reference to the original source (see Ref. [6 ). The same disadvantage char- 
acterizes the results of other papers [2HI EHj because the exact solution [HI] is not familiar 
to the authors [95ll9fij . 

Here we shall use the functional formulation of NBG (see, e.g., Ref. |3D] ) which seems 
to be more convenient for RG studies [7^1 !H]- We shall be mainly interested in the 
zero-temperature critical behaviour of NBG [HIj. The relation to other model systems 
with similar properties and recent research will be also mentioned (Sec. 5.7). The func- 
tional NBG formulation is given by Eqs. (27) and (28) where the parameter r = — \i is 
related to the (bare) unrenormalized chemical potential /i < 0. In order to avoid lengthy 
mathematical formulae we shall set o = 2. The generalization of results to < a < 2 
is not difficult. The critical behaviour is investigated with respect to variations of the 
chemical potential \i from the value /i c = /i(T c ), corresponding to the critical point (for 
IBG, /i c = 0). A treatment of the phase transition in NBG at a constant density with the 
help of variations of the quantity t ~ (T — T c ) has been performed in Ref. |85j . 

5.2. Renormalization group equations 



The RG recursion relations within the one-loop approximation (to first order in e = du—d) 
are [8T] : 



K = bk, 




m 



(38) 



r 



r + -(n + 2)(v/P)I 1 (r) 



(39) 



and 



(y/P)> = b*-^Uv//3)-±(v/P) 2 [(n + 6Mr) + 2I 2 (r) 



(40) 
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Here we shall write in an explicit form the quantities 



h{r) =(3K d J dkk d - l n{k), (41) 



A/6 



A 

/A r , f „ k d - 1 



Ur) = £ HQ / dfc^- — — — , (42) 

A/6 

and 

/ 2W = §^/\^-«^. (43) 

A/6 

In the above expressions the integration is performed up to an upper cutoff A <C (vr/a); 
b > 1 is the RG rescaling factor [TJ H2] . 

In contrast to other RG studies here the problem for the cutoff A is not trivial [23 • 
The problem arises in the treatment of the dilute (d > 2) -dimensional Bose gas in the 
low temperature limit (T c — ► 0) which corresponds to an extreme dilution p —>■ 0; see 
Eq. (23) which indicates that A ~ (tt/o) — > for p = a~ d — ► 0. It seems at first sight 
that the interparticle distance a cannot be used in the definition of the finite cutoff A 
or as a length unit of dimensionless characteristic lengths £/a and A/a. For this reason, 
the fourth (and last) characteristic length in the problem - the finite scattering length 
v ~ {h 2 /m)a sc iniHD] was discussed as a possible inverse cutoff, (1/A) ~ a sc , for the first 
time in Ref. jHI] and later used also in Ref. |2S1; for d ^ 3, v ~ (h 2 /m)aii 2 . 

Let us note, however, that the RG investigation can be reliably performed by the standard 
cutoff A = n/a, although a tends to infinity for the zero temperature critical behaviour. 
In fact, according to Eq. (23), A > a for all temperatures in the interval < T < 
A((d/2)T c (p) including the close vicinity of T c where £ > a, too; note that A((d/2) > 1. 
In the temperature range of interest the characteristic lengths £ and A are always greater 
than a and, therefore, there is no danger of shortcomings in the description; see also the 
brief discussion after Eq. (50) in Sec. 5.3. 

The summation over the frequencies ui of the internal lines of perturbation diagrams 
formed by the legs of one and the same Hamiltonian is performed with the help of the 
rule 

J>(+°)^G (g) , (44) 

where (+0) denotes a positive number which is set equal to zero after the calculation. 
The correlation function Go(q) is given by Eq. (28); c.f. Refs. PUEJIBB!, where the Green 
function Go(q) has an opposite sign. 
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The summation (44) leads to a difference between the integrand n(k) in Ji(r) and the 
corresponding integrand in the RG relations presented in Refs. [HE] but coincides with the 
corresponding result in Refs. jHJEElIB]- Moreover, the integrand in I^ir) differs from 
the corresponding integrand in the RG equations in Ref. [7H [75] but coincides with that 
in Ref. [HE] as mentioned for the first time in Ref. [SI]. The mentioned errors in papers 
before 1981 do not have an effect on the results for the critical behaviour in the classical 
limit (g < 1,T C ^ 0). However, in the low-temperature limiting case (T c — > 0), these 
errors lead to wrong results for the FP coordinates and the rational extension of the loop 
expansion in higher than first order. 

Owing to this formulation of the RG treatment the relation (38) for ui should be referred 
to the temperature: 

T' = b 2 - v T . (45) 

For the same reason the mass m should be kept invariant. We shall proceed with this 
choice up to the end of this subsection, where we shall propose a more elegant scheme 
of scaling. The mass (m— ) invariance implies 77 = 0. This is the usual result within the 
one-loop approximation. Apart from quite special cases the exponent value r\ = receives 
e-corrections of order 0(e 2 ), i. e., in the two- and higher-loop approximations. 

The task is to solve RG Eqs. (38)-(40). This means to obtain FPs and investigate their 
stability. Eq. (45) shows two temperature FP values: Tq = 00 (classical) and Tq = 
(quantum). The infinite FP value Tq of the temperature has nothing in common 
with infinite temperatures. Rather it describes the classical critical behaviour near finite 
temperature critical points T c > 0. 

Remember that the classical limit of quantum Hamiltonians (27) is taken, in a strict 
mathematical sense, by decreasing the interval [0, /3] of variations of Matsubara time r to 
zero which corresponds to the limit T — > 00. Obviously, the FP coordinate Tq = should 
correspond to the low-temperature regime. As FPs lie on the critical surface [£(T C ) = 00] 
in the Hamiltonian parameter space it is clear that the zero FP value corresponds to 
a zero-temperature critical point (T c = 0). However, it can be shown that for "high- 
temperature" critical points the parameter T is absolutely redundant together with the 
relation (45) and the corresponding infinite FP value of temperature (see Sec. 5.5). The 
infinite FP value T£ = 00 merely indicates that the zero- temperature critical behaviour 
is unstable towards temperature fluctuations. 

Eq. (38)-(40) can be solved analytically in two limiting cases: 

(i) (3e{k) 1 (high-temperature behaviour), 

(ii) f3e(k) ^> 1 (low-temperature behaviour). 

These conditions will be discussed in Sec. 5.4. Here we shall take the leading terms of the 
integrals (41)-(43) for each of these limiting cases and perform the formal analysis of RG 
equations. 
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The limit (i) yields the classical universality as shown in Refs. |7l^ 175] . In this case 
(e = 4 — d) the FP value of the temperature is Tq = oo. This allows small values of the 
factor p in the condition (i). The RG relations have two stable FPs. For dimensionalities 
d > 4, the Gaussian FP with (r, u)-coordinates r G = w G = is stable and describes the 
usual MF behaviour; u = v/(3. For dimensionalities 2 < d < 4, the Heisenberg FP is 
stable. The coordinates of HFP are given by 

(n + 2) h 2 A 2 16tt 2 /ft 2 V 

TH = ~WT^)^ UH = ^Ts) l^J e • (46) 

As the temperature FP coordinate is infinity the FP value vjj of the parameter v = /3u is 
equal to zero because the respective FP value % of u is finite, as given by Eq. (46). This 
behaviour of the interaction parameter (v or u) is not strange because it precisely describes 
the behaviour of the system at high temperatures where the modes 0(g) with u>i ^ can 
be neglected. Ignoring the quantum fluctuations we obtain a classical Hamiltonian with 
an interaction constant of the form u = v / (3 which plays the role of an interaction constant 
in the usual classical Hamiltonians (7i = f3H). 

However, this is not the final answer. As shown in Sec. 4, we must substitute the modes 
0(0, /c) with the true classical modes \/]3<f>(k) and this yields a factor (3 ~ (l/k B T c ) in 
front of parameters in the 2 -part as well as the square of the same factor in front of 
the interaction constant. So, in the classical variant of the theory the actual interaction 
constant is (3v. Therefore, the parameter u = (v/f3) is the interaction in the classical 
Hamiltonian in terms of the field 0(0, k) and this parameter transforms to u = j3v when 
the classical Hamiltonian is written by the field <p(k) = 0(0, k)/y/]3. 

This discussion can be used to explain the infinite value of FP coordinates of temperature: 
Tq = T^j = Tq = oo. The field 0(0, k) has a rescaling factor b (for rj = 0). The field 
crossover [0(0, k) = -^P(j)(k)] implies that this rescaling factor should be associated with 
the proper classical field 4>(k) rather than with the temperature. This is just what happens 
in the classical limit (i) where the relation (45) does not exist at all. Therefore, in the 
limiting case (i) the temperature is a redundant parameter and the relation (45) can be 
suspended from the high-temperature analysis although it remains in the general RG 
scheme. Up to the stage, for which the investigation is performed in terms of the field 
0(0, k), i.e., in the low-temperature (classical or quantum) regime we must consider the 
relation (45) and the FP temperature value Tq = 0. 

5.3. Low temperature behaviour 

In the low-temperature limiting case (ii), Eqs. (38)-(40) are solved with the help of e = 
(2 — (i)-expansion. This expansion reflects the dimensional CQC which is given by D = 
(d + 2). The upper critical dimensionality dy is changed from du = 4 for the classical case 

(i) to du = 2 for the quantum case (ii) |20j. Up to now there is no evidence that the case 

(ii) describes quantum critical phenomena but it is clear that the low-temperature limit 
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can be taken within the RG scheme and that it will bring to a new critical behaviour. 
Within RG this limit always exists because of the lower cutoff (A/6) > which permits 
the condition (ii) irrespective of the value of the ratio (r/ksT); see also Eq. (21). For 
the case (ii) the integral -?2( r ) tends exponentially to zero but the integral ^{r) is finite 
and exhibits a power-law infrared divergence for d < 2. This yields the borderline value 
djj = 2. Performing standard calculations we obtain Eqs. (39) and (40) in a simple form 

r' = b 2 r , v' = b 2 ~ d v(l - a v) , (47) 

where 



a = 



h(r) 



A 

1 f k dk 



Pe(k)^oo An J e(k) 

A/b 



(4£ 



The straightforward calculation gives 



a = — Inb . 49 

The relevant parameters are T and r. The variations of these parameters near the FP 
values T* = r* = drive the system away from the zero-temperature critical state. The 
inclusion of the temperature T as a second relevant parameter corresponds to a real 
physical situation, namely, that the zero-temperature critical state is approached when 
both T and r tend to zero. The zero- temperature critical behaviour is described by two 
stable FPs: the Gaussian FP (Tq = tq = vq = 0), which is stable for d > 2, and the 
Gaussian-like FP [BT| . 

t n 27Th2 

T G \ = r G \ = , vqi = e . (50) 

m 

The equation for vqi shows that the renormalized value of s-wave scattering length is 
a sc ~ (l/e) 1 ^; e = (2 — d). The scattering length a sc — > oo for the important case e — > 
of 2d Bose fluids and, therefore, the cutoff A ~ (l/a sc ) applied in Ref. |2S1 seems to be 
inconvenient. 

In the remainder of this Section we shall discuss the nontrivial Gaussian-like FP which 
will be called G1FP. GIFP has Gaussian values for the main critical exponents: rj = 0, 
v — 1/2, z — 2. But the interaction parameter gives a correction-to-scaling exponent and 
this circumstance does not allow to put this critical behaviour in the Gaussian universality 
class (such an incorrect conclusion has been made in Ref. jHZj). Note, that the Gaussian 
universality class is sometimes denoted by n = — 2 (see, e.g., Ref. |72*j). 

In the present case the e-analysis can be extended to any order in e [81 ^ and this excep- 
tional case will be discussed below. The problem has an exact solution jHIj because of the 
great simplification of the perturbation series in the limit (ii). In this limit, all self-energy 
perturbation contributions (see Fig. 3b) exponentially tend to zero. In each diagrammatic 
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Figure 3: (a) A tree diagram denoting the interaction part of the Bose Hamiltonian. (b) 
The compact self-energy diagram which is equal to zero in the limit T — > [the thick 
loop denotes the full (renormalized) Green function G(q)]. (c) An example of a diagram 
from the perturbation series for the interaction vertex which gives a zero contribution in 
the zero temperature limit, (d) The infinite ladder series of diagrams which yields the 
geometric progression (51). 

(perturbation) term in the self-energy function (TJ ED] all poles in the frequency (a>) 
lie in the complex upper-half plane and, hence, the contour of the respective frequency 
integral can be deformed into the lower-half plane to give zero. Therefore, in this limiting 
case, no corrections to the MF values of the critical exponents describing the main scaling 
laws will appear to any order of the loop expansion. 

Furthermore, in the same zero-temperature limit (ii) the bigger part of perturbation 
contributions to the interaction vertex v' tend to zero, too. For example, the diagrams 
of type shown in Fig.3c do not give contribution. The only exception is the ladder series 
shown in Fig. 3d. This is the so-called superconductivity channel for the interaction 
vertex where the Green function lines are oriented in one and the same direction; for 
more details, see Refs. [TJ |HH] • The final result of the summation of the ladder in Fig. 3d 
is an infinite geometric progression [81J 

1 + a v 

The above consideration proves that within the grand canonical formalism the dynamical 
and static critical exponents of main scaling behaviour described by G1FP have Gaussian 
values. This is true for the scaling laws written in terms of the chemical potential \fi\ = r. 
As the parameter r does not receive a renormalization it remains equal to that of IBG. 
The zero temperature critical point is defined by [r(T c )/T c ] = as is for IBG (Sec. 3). 

The FP value Vq\ can be obtained from Eq. (51) to any order in e = (2 — d). To do this one 
should expand the integral (48) and Eq. (51) to the corresponding order in e by having 
in mind that the FP value i>gi(c) have to be also expanded. For the sake of convenience 
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in calculations to second order in e, the FP value VQ\(e) is usually written in the form 



VGi(e) = ^-[1 + elnA] (52) 



Expanding K d in e we have 



2tt 2 . 

v gM = — M 1 + e 



'^ + ln- A 



2 v 7 ^. 



(53) 



(ce is the Euler constant). The FP coordinate vgi receives corrections in all orders of the 
e expansion. Within the one-loop approximation, the stability exponent y v associated 
with the interaction parameter v has a negative value y v = — e = (d — 2) for < d < 2. 
This exponent does not receive higher order e corrections and, therefore, the value (d — 2) 
is exact. 

This completes the proof [HIj. The above results have been confirmed by Monte Carlo 
calculations for one-dimensional interacting bosons |97j . 

We should emphasize that G1FP is non-Gaussian FP. It is a product of the thermal and 
quantum fluctuation interactions represented by the interaction parameter v. In fact, 
G1FP is conjugate to actual GFP (vq = 0) within the same RG analysis. Moreover, the 
total critical behaviour includes corrections to the main scaling laws which are described 
by G1FP. The interpretation jHZ| of results within the classical universality class (d + 
2, n = —2) has no heuristic significance. As we have already mentioned this interpretation 
is incorrect because of the correction- to-scaling exponent coming from the interaction. 
Therefore, G1FP constitutes a new universality class of zero-temperature QCP which has 
numerous applications, in particular, to fluids of real bosonic atoms and XY magnets 
(Sec. 5.7). 

5.4- Breakdown of quantum universality 

G1FP and the critical behaviour described by it are quite unusual. Let us mention that 
here the term "usual zero-temperature FP" means a zero-temperature FP which can 
be obtained from the respective stable finite temperature FP of the same system after 
substituting e = (djj — d) with e = (djj — d — z ) in the formulae for the coordinates and the 
attached critical exponents of the latter. Therefore, there is a clear form of conformity 
between the finite temperature FP and the respective zero-temperature FP which we 
temporally call "usual" zero-temperature FP. This point of view is consistent with the 
definition of universality of quantum critical phenomena given in Sec. 4.3. 

If all zero-temperature FP are usual in the above mentioned sense then one does not 
need to define new universality classes for quantum critical phenomena. The universality 
classes determined by the known finite-temperature FPs can be used, after the account 
of the dimensional shift due to the CQC, to describe the possible types of quantum 
critical phenomena, too. If this scheme works, we can introduce the term "universality" 
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of quantum critical phenomena by understanding this as a quantum universality which 
coincides with the classical universality at a shifted effective dimensionality of the space 
(see also Sec. 4.3). Then a zero-temperature FP obeying the mentioned correspondence 
with the respective finite-temperature FP can be called "FP describing quantum critical 
universality" or, shortly, "quantum universality FP." 

Certainly, the G1FP is not such usual or quantum universality FP. It does not obey 
the present definition of quantum universality although it comes as a result of CQC in 
NBG. According to Fisher and Hohenberg [HH] the G1FP describes a "quasiuniversality" 
associated with the existence of an almost marginal line in the Hamiltonian parameter 
space (T, fi, v) defined by T = fi = and a nonzero interaction parameter v which can 
take positive values. 

This is in contrast with other zero-temperature FP revealed by J. Hertz j20j which are 
usual in the above mentioned aspect and obey the quantum universality. Therefore, there 
is a number of real systems which do not satisfy the widely accepted notion that the 
main effect of the quantum fluctuations on the zero-temperature critical behaviour is 
CQC. Such systems are the Bose fluids and the XY magnets (Sec. 7) where the quantum 
universality concept is not valid. 

5.5. Discussion of limiting cases 

The limiting cases (i) and (ii) are analyzed under the general condition (£/A) ^> 1. For 
a microscopic model, A ~ (tt/o) and, hence, we have (£/a) 1. Then the inequality (i) 
yields A < a, i.e., the high-temperature condition. The values A < d have a statistical 
meaning in the continuum limit (V/a d ) — > oo. 

The high-temperature condition is certainly satisfied near the infinite FP value Tq = oo 
and well below it. The condition (ii) is valid either in the classical low-temperature critical 
region 

a < A < f , (54) 

or in quantum one, given by 

a < f < A ; (55) 

c.f., the general criterion (3). The inequalities (54) and (55) can be written in terms of 
parameters T and r which enter in the RG relations. The condition (54) will be valid 
near and at low-temperature FPs given by the zero values of T and r provided in the 
zero-temperature limit T — > we have (r/ksT) — > 0. In the same limit, the quantum 
condition (55) will be valid, if (r/ksT) — > oo. In particular, to answer the question which 
of these two types of the low-temperature critical behaviour is described by G1FP (50), we 
must know the value of the ratio (tqi/Tqi). If it is greater than unity, G1FP will describe 
quantum critical phenomena but if it is less than unity, the low dimensional (d < 2) 
low-temperature critical behaviour will be classical. 

A standard example is given by the system of interacting bosons at a constant density (the 
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case of noninteracting particles has been discussed in Sec. 3.4). As the perturbation series 
does not yield self-energy contributions at all the parameter r for ideal and interacting 
Bose systems is the same. The value of the correlation length critical exponent u, cited in 
Table 1, is greater than the exponent 1/a of thermal length A. So in the zero-temperature 
limit the ratio (r/A^T) ~ (A/£) 2 will tend to zero; c.f. Eq. (21). Therefore, the low- 
temperature critical behaviour described by G1FP will be classical which is true for T- 
driven transitions (Sec. 2.3). The p-driven transitions will exhibit quantum critical effects. 
These results are consistent with the general criterion (9). 

In the same way one can show that the critical effects near the phase transition points of 
interacting bosons under the condition of a very low constant pressure will be influenced 
by quantum effects (see also Sec. 3. 3). The application of the results from Sec. 5.3 to dilute 
Bose systems in the low density limit p — > will be briefly discussed in Sec. 5.7. 

The NBG properties at finite low temperatures are very similar to those of IBG at constant 
density (Sec. 3.4). It is then convenient to investigate the corrections to the (spherical) 
Hartree limit by the standard 1/n— expansion J2|- The alternative approach j37j by a 
reduced Landau-Ginzburg free energy function has been mentioned in Sec. 4.1. 

5.6. Formulation by thermal wavelength 

An alternative RG treatment can be performed by the transformation 

^H 2 /2m0(q) -> <p(q) (56) 

in the Bose Hamiltonian (27)-(28). In the new notations for the theory parameters 
the mass m is absent and the Matsubara frequency ui is substituted by 87r 2 Z/A 2 , i.e., 
by (1/A 2 ). The parameters r and v are multiplied by factors (2m/h 2 ) and (2m/h 2 ) 2 , 
respectively. The important point is that the recursion relation (38) for the mass m and 
the recursion relation (45) for the temperature are now substituted by 

1 = b\ (57) 

and 

A' = r x A , (58) 
respectively. Eq. (58) can be written in the form 

m'T = b 2 mT . (59) 

The comparison of Eq. (38) for m and Eq. (45) for T with Eqs. (57) and (59) shows that 
the second formulation is not completely equivalent to the usual one although the formal 
results for the critical behaviour are equivalent. Within the present formulation Eqs. (57) 
and (59) do not imply the mass invariance. Rather the mass enters as a factor in the 
parameters A, r and v and, hence, participates in the renormalization scheme but in a way 
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different from the suggested in Ref. [71^ 175] . From Eq. (59) one cannot conclude anything 
about the individual behaviour of T and m towards the RG transformation. Rather the 
correct conclusion is that the thermal wavelength A is renormalized, as given by Eq. (58), 
and that RG transformation will drive A to zero, unless it is at the high-temperature FP 
value Xq = 0; another FP value of A is \q = oo. 

Note that the renormalized thermal wavelength A = A (b) does not belong to the para- 
metric space (r, v ) of the Hamiltonian. Hence, the temperature T is also excluded from 
this space. It is convenient to consider these parameters as describing HLTC. The renor- 
malized wavelength A (6) drives the system from the high temperatures (A <C a) to the 
zero temperature (A 3> a) Gaussian-like critical behaviour; see Eq. (58). The rescaling 
factor cannot be taken as a crossover parameter because of the RG restriction e ln6 < 1. 
So the RG flows from the usual Heisenberg FP where AJ. = to G1FP where Xq = oo 
are produced by successive RG transformations. If the temperature is equal to zero, the 
system will be exactly at G1FP and the FP Hamiltonian is invariant towards RG. The 
system remains in the ground state of total BEC (po = p)- For any T > 0, the RG flow 
will drive the system to the usual classical behaviour. Now it is not difficult to rederive 
all RG results in the present formulation by avoiding the unnecessary mass-invariance 
condition and the unnatural HT FP value Tq = oo. 

5.7. Related topics and applications 

As the systems undergoing zero-temperature phase transitions are numerous here we shall 
mainly discuss the applicability of the results from Sec. 5.3 to real systems and topics 
having some relationship with these results. 

5.7.1. Magnetic and ferroelectric systems 

Experiments jHHl HUH H00| HOlj on quasi-one- dimensional antiferromagnets with integer 
spin in an external magnetic field H, mainly intended to a search of Haldane gap |102j , 
reveal a continuous phase transition to a nonzero ground state magnetization which oc- 
curs at the absolute zero (T c = 0,H C > 0). Theoretical investigations |103| 11041 H05j 
demonstrate that this phenomenon is a type of BEC of magnons with spin 5 = 1. The 
critical properties in the vicinity of the transition point (0, H c ) are described by G1FP 
(Sec. 5.3), as shown in Ref. [Hi], where low-temperature (T ~ 0) features and the phase 
diagram of these spin chain antiferromagnets have been investigated. 

The results for interacting real bosons are straightforwardly extended to systems described 
by the quantum XY model in a transverse field and, in particular, to XY magnets without 
a time reversal invariance. RG calculations to the one-loop order in Refs. [7T| IT2] have 
drawn the attention to the possibility for a Gaussian zero-temperature critical behaviour 
of this model. In contrast to systems with different symmetry, for example, uniaxial 
ferromagnets and ferroelectrics described the TIM (Sec. 6), in XY systems above the 
critical transverse field the spins at T = have maximal projections on the field axis and 
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there are no quantum spin fluctiations at all. Below the critical field, in the ferromagnetic 
phase, such quantum fluctuations exist; for the ground state properties of spin systems 
see Ref. jHHj. 

In Ref. [71] the real symmetry of XY systems has been neglected and the obtained RG 
equations depart from the correct Eqs. (38) - (43) just as the RG equations in Ref. [Ho] (see 
Sec. 5.2). The ignored symmetry properties are related to the particular form of the bare 
correlation function Go(q) which has been obtained in Refs. [72~ 1 11071 llt)Sj exactly of the 
form (28). The main interaction term in the effective Hamiltonian of the XY model is of 
the form given by Eq. (27), but there are also additional interaction terms of type |0| 2m ; 
m > 2. It has been shown in Ref. |109j . that the additional interaction terms do not 
change the critical behaviour predicted with the help of the 4 -Hamiltonian (27)-(28) 
and, therefore, the RG results from Sees. 5.2-5.6, can be straightforwardly applied to XY 
systems; see, e.g., Ref. [73~] . 

We should emphasize that the transverse field in the XY model plays the role of the auxil- 
iary parameter X, introduced in Sec. 2. The X— driven transitions considered in Sec. 2.4 
are easily performed in XY systems by variations of the transverse field around its critical 
value at a fixed low temperature. Similar field variations in uniaxial ferroelectrics [T8* llll0j 
and ferro magnets PU Ej described by TIM are discussed in Sec. 6. The measurements 
in real substances should indicate the Gaussain-like critical behaviour in XY systems 
and the standard zero-temperature critical behaviour described by the universality Ising 
quantum FP of TIM [20] • The results for the quantum ctitical behaviour described by 
GLFP have been applied [HI] also to low- dimensional quantum antiferromagnets. 

Therefore, according to our considerations in Sec. 2 as well as the concrete studies of 
nonuniversal characteristics of XY systems |1061 11071 HH| - quantum critical phenomena 
should be observed in low-temperature critical experiments on low dimensional XY fer- 
romagnets. The arguments presented so far demonstrate that the low-temperature phase 
transitions in XY systems are a promising area of experimental studies of the Gaussian- 
like critical behaviour. Besides, respective experiments on uniaxial ferromagnets and 
ferroelectrics described by TIM are expected to confirm the universal quantum critical 
behaviour. 

5.7.2. Quantum Hall liquid and superconductivity 

The results from Sec. 5.3 can be applied to the (fractional) quantum Hall liquids (see. 
e.g., f52ll!13j . This problem has been investigated by Schakel |114|lll5j : see, also Ref. jlj. 
The quantum Hall liquid was considered |1141 I115j in the framework of the effective 
(quasimacroscopic) Chern-Simons-Ginzburg-Landau (CSGL) theory which describes the 
fractional quantum Hall effect (FQHE) [221 H13j . This effective field theory is based on 
the 4 — theory given by Eqs. (27) and (28) but the difference is that the Bose field <p(x) 
interacts with the sum of the vector potential A(r) of the external magnetic field H and 
the Chern-Simons gauge field which is related to internal degrees of freedom (composite 



41 



particles featuring the FQHE). It has been shown that the basic features of the quantum 
phase transition in this system are described with the help of the exact solution presented 
in Sec. 5.3 (see also Sec. 7.6.2 for a discussion of disorder effects). 

The relationship between the superconductivity and the quantum criticality seems to 
be a new important issue in strongly correlated materials, including heavy-fermion and 
high— T c cuprate superconductors, in particular, in cases of unconventional d— and p— wave 
Cooper pairs (see, e.g., Refs. HIT] ). The problem has been mentioned for the first 
time within the framework of the resonating-valence-bond theory of high-temperature su- 
perconductors |118| 1119] ; see also Ref . for a discussion of quantum critical points in 
superconductors. Here we shall draw the attention to several new results. 

Barzyukin and Gor'kov |12()j have recently shown that in two-dimensional superconduc- 
tors, the phase transition line to the Larkin-Ovchinnikov-Fulde-Ferrell state may extend 
up to T c = at a certain value of the external magnetic field. If the theory is reliable, 
this phase transition line may give a new type of quantum (multi) critical point, or, a 
zero-temperature first-order transition. 

Other magnetic field induced quantum phase transitions were established experimen- 
tally in heavy-fermion superconductors such as Ce-based (CeCu6-xAg x |121| I122j and 
CeCoIn 5 |123j ) and YbRh 2 Si2 |124j compounds (see also Ref. |125j ). The bilayer ruthen- 
ate Sr 3 Ru207 exhibits a low temperature metamagnetism 126 . Assuming that for a 
certain value of the external magnetic field the phase transition line may reach the ab- 
solute zero, theoretical arguments |127j have been presented in favour of a weakly first 
order quantum phase transiton in Sr 3 Ru207. 

A phenomenological model |128j . that proposes a coexistence of the ferromagnetism and 
the unconventional superconducting state of spin-triplet Cooper pairs in heavy-fermion 
compounds, yields, in accord with experiments |129| I13U| 1131 j . a two pressure-driven 
quantum phase transitions at the same zero-temperature multicritical point: a transi- 
tion from ferromagnetic to paramagnetic phase and a transition from superconducting to 
normal metal. 

The theory j!32j and the experiment |133j on ultrathin superconducting wires demonstrate 
the crossover from thermally activated phase slips (a product of thermal fluctuations) 
and quantum phase slips, activated by the quantum fluctuations. This crossover very 
much resembles the general scaling picture of effective HLTC and, in particular, CQC at 
nonzero temperature, outside an unattainable by experiment, small asymptotic vicinity 
of the critical point (see the discussion in Sec. 2). Dissipation-driven quantum phase 
transitons in (quasi-) one- dimensional arrays of Josephson junctions |134j are a subject of 
intensive experimental [13BJ and theoretical |136[ I137j research. 

A Monte Carlo calculation [13.8! based on a two-dimensional boson Hubbard model, equiv- 
alent to an anisotropic spin-1/2 XXZ model with definite parameters, reveals a phase dia- 
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gram with a quantum (T = 0) superfluid-solid phase transition. As mentioned in Ref. j!38j 
there is a striking qualitative similarity of this phase diagram to those of fermionic 3 He 
and bosonic 4 He two-dimensional systems. 

5.7.3. Dilute Bose fluids 

At the end of this discussion let us come back to the systems of real boson particles. 
A potential application of the nonuniversal critical behaviour described in Sec. 5.3 is 
to thin films of 4 He |139[ I14()j . An interesting opportunity exists for an experimental 
investigation of the crossover of the critical behaviour of IBG to the critical behaviour of 
NBG by a careful controlled overall density p of bosons up to the dilute limit (p <C 1) in 
4 He g!41| 11421 H43t H441 1145] and spin-polarized hydrogen j!46j . The critical fluctuations 
of the superfluid order parameter < <f>(x) > can be accounted by the experimental data 
for the behaviour of the superfluid density p s ~< |0 2 | >. 

Experiments following this idea have been performed |141^ 1142] for 4 He condensed in Vicor 
which is a highly connected semiregular sponge-like glass. In other experiments |143| 11441 
I145j another porous medium (Silica Gel) has been used. This is again a medium which 
consists of randomly distributed three-dimensional (d = 3) networks of interconnected 
pores and channels of diameter approximately 0.005 pm (the pore size can be varied from 
sample to sample). In these experiments the critical temperature T c of the superfluid 
phase transition is observed to decrease with the decrease of the density p of 4 He. The 
respective curves can be reliably extrapolated to zero T c (p = p c ) for a nonzero critical 
density (p c > 0). 

The experiments |141[ I142| I143| 11441 I145j have been extensively analyzed in theoretical 
works [23 121 HIl Q3E1 QUI EUl EH based on scaling and RG methods. All mentioned 
porous media are a source of disorder, so that certain properties of the Bose fluid in such 
media should be described with the methods of the theory of disordered systems (Sec. 7). 
However, some main features can be discussed by the means of pure systems where no 
disorder effects are present. 

As pointed out in Ref. [21] such systems undergo a superfluid onset transition as the 
density p increases at T = and, moreover, this phase transition may be related with 
the theory presented in the preceding Sees. 5.2-5.6. A similar conclusion has been made 
in Refs. [11481 EH] where the empirical data is interpreted with the help of the crossover 
from the critical behaviour of NBG to that of IBG. These notes demonstrate that the 
experiments are close to a confirmation of the Gaussian-like critical behaviour established 
in Sec. 5.3, provided the respective porous medium can be considered as a factor leading 
to a decrease of the (effective) spatial dimensionality of the Bose fluid contained in it up 
to d < 2. Similar experiments can be made on thin superfluid films. 

Following Ref. [21] we shall briefly mention some properties of this experimentally ob- 
served phase transition in 4 He placed in Vycor. While the experimental feasibility of 
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the probing low-temperature behaviour stops short of T = 0, the superfluid transition at 
sufficiently low temperatures will be dominated by GLFP until the very close vicinity of 
the critical density p c is approached. The crossover from the zero-temperature critical be- 
haviour described by G1FP and the finite-temperature critical behaviour corresponding to 
T c > 0, described by the standard XY FP(n = 2) has been treated in Ref. [21] as a stan- 
dard multicritical-to-critical crossover (within this interpretation the zero-temperature 
critical point is considered as a multicritical point). 

We should draw the attention to a quantum (T = 0) phase transition |15f j to superfluid 
state in 3 He in aerogel at a nonzero pressure and a critical density of 3 He p c > 0. Unlike 
bulk 3 He which is superfluid at all pressures (densities) between zero and the melting 
pressure, 3 He in aerogel will be not superfluid, unless the density of 3 He exceeds a critical 
value p c > 0. Some of the problems of the quantum phase transition revealed in Ref. |15T| 
are similar to that for the respective phase transition in 4 He, discussed in a few lines above 
(for details, see Ref. |151j ). 

6. Transverse Ising model 

6.1. Microscopic and field models 

Here we shall consider the quantum transverse Ising model with the aim to reveal the 
properties of HLTC. We shall use MF, the lowest-order perturbation theory and Ginzburg 
criterion considerations and RG. 

TIM describes essential features of phase transitions in ferromagnets with a strong uni- 
axial anisotropy |152j and displacive phase transitions in certain types of quantum fer- 
roelectrics [TH1 1153| I154| I155| IllOj : for experiments, see, e.g., Refs. j2"""[ I156j (for other 
applications of TIM, see Refs. [1571 I15£"| ). Here we shall mention the MF results |152j . 
the cumulant expansion |157j . the RG investigations [T""". |2"U1 ITo"4*| I155j . the description of 
CQC in Hartree limit [7UJ 1159j and the formal interrelationship between the finite-size 
crossover in systems of slab geometry and CQC in TIM [70J. 

TIM is given by the Hamiltonian j!52j 



where 5 7 , 7 = (x,y,z), are the components of spins with a magnitude S = 1, is the 
exchange interaction and T is the transverse field magnitude. We shall assume that TIM 
is defined on a <i-dimensional regular lattice with a lattice spacing a = 1. The assumption 
for nearest-neighbour (nn) interactions which we shall use to define the parameters of 
the effective Hamiltonian, does not restrict the generality of results in this Sec. 

The fluctuation field Hamiltonian of TIM can be derived from the microscopic spin 
model (60) by the Hubbard- Stratonovich transformation [T""" [ 170* 1 This Hamiltonian 




(60) 
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(7i = H/T; ks = 1) is given in the form 



a,c/ 



where 



«,/3;<?i<?2<?3 

is a real scalar field with a (bare) correlation function in the form 
G \q) = \u l \ 2 + k 2 + t . 



The parameters to an d «o are given by the expressions, 
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In Eq. (64), J = J(0) = 2<iJo is a product of the number z = 2d of nearest neighbour nn 
spins and the constant Jo of the single (i — j) nn interaction (Jo = J%j for nn sites i and 

It is convenient to use units, in which all quantities in the effective Hamiltonian Ti, are 
dimensionless. The dimensionless wave components ki are given by fcj = 2'KKl i /L il where 



Jth(r/T) 



2dT 



(65) 



In the critical region to < 1, k ~ 1. The upper cutoff for the wave numbers k is A = 77r/t, 
where 7 is a small number (7 <C 1) |4*5] . 

The characteristic lengths ^ and A are given by £ = 1a/|£o| and A = (T/T). We shall 
consider only the low temperature domain defined by A ^> a. In this case, u = (J 2 T/8T 3 ). 

6.2. Mean field approximation 

Within the standard MF approximation we shall consider only the uniform mode 0(0) of 
the field 0, that is the q— dependent (classical and quantum) fluctuations will be ignored. 
It will be more convenient to use the MF free energy Q = T7^[0(O)] instead of the 
dimensionless MF free energy 7i[0(O)]. This choice, after a change of the nonequilibrium 
order parameter from 0(0) to 0o = (T/V) 1 ' ' 2 0(0), will induce an extra factor 1/T in 
front of 0Q-term of free energy Q which makes possible to avoid difficulties in our further 
analysis connected with the definition of the order parameter 0o at T = 0. The problem 
is that the parameter m is proportional to T, see Eq. (64), and the investigation with 
the original parameter 0(0) will give for the equilibrium order parameter 0(0) ~ (1/T) 1 / 2 
which is divergent for T — > 0. This is an example of order parameter HLTC 



The Gibbs thermodynamic potential in the form 

Q = V 



|*8 + f*S 



(66) 
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Figure 4: (a) The graphical representation of Eq. (67) for J = 1. (b) The high- and 
low-temperature parts of the curve in Fig.4a with an indication of T— and T— transitions 
(see the text). 

allows a correct MF analysis at low temperatures. The analysis of the free energy (66) 
can be done straightforwardly and we shall not enter in details. The critical line (Fig. 4a) 
is defined by t (T c , T) = 0, i.e., 

T << r) = Amu) (67) 

or, equivalently, by t (T,T c ) = 0, which yields 

£ = «*> 

Despite of the simple form of Eq. (66) the MF critical properties of TIM cannot be inves- 
tigated analytically for the whole curve T C (T) because of the quite complex dependence 
of parameters t and u on T and T. 

We shall consider the low temperature critical behaviour in the vicinity (|T — J| < J) of 
zero temperature critical point [T c ( J) = 0, T c (0) = J] and the neighbour critical points 
with coordinates T c J and T c ~ J. In this low temperature region we can distinguish 
four types of phase transitions along the lines: 0J, AJ, aa, and bb; see Fig. 4b. We 
shall suppose that the couples of parallel lines are very near to each other. Following the 
terms introduced in Sec. 2, the transitions along the lines A J and aa are T-transitions 
whereas the transitions along the lines 0J and bb can be thought of as T-transitions. The 
results can be compared with those for the high temperature T— and X— transitions in 
the vicinity of the critical point T c (0); see the respective lines 0J, bb, AJ, and aa in Fig.4b. 

In MF we shall consider the critical exponents j3, 7, and is, of the order parameter |0o| = 
(Tltol/4-uo) 1 ^ 2 , the susceptibility x — VI*o|, and the correlation length £ = 1 / 1 1 1//2 , 
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respectively, where t = t for t > 0, and t — 2to f° r t < 0. Obviously, the behaviour of 
these quantities with respect to variations of T and T depends on the behaviour of to- Note, 
that parameter to as given by Eq. (63) does not change in the Gaussian approximation for 
the quantum and classical fluctuations. This makes possible to extend the consideration 
of the parameter t in the paraphase (0 O — 0), where the results should be interpreted as 
pure fluctuation effects. 

The summary of results reads: 

Line 0J: We obtain t ~ (r — J)/T, O ~ l^o| 1//2 > X ~ r/|T — J| and, therefore, 7 = 1 
and f3 = 1/2. There exists a full coincidence between the standard MF behaviour with 
respect to variations of T around J at T = and the present critical behaviour due to 
variations of T around J at T = 0. This correspondence can be written as T <-> T. 



Line bb: The parameter to is 



^^(1 - ^e- 2r * /T + ..-)■ (69) 



Here to tends to the value (r — J)/Y when T C (T) — > J for T — > 0. The critical exponents 
are the same as those along the line 0J. The exponential correction in Eq. (69) can be 
neglected. 

Line A J: There is no equilibrium ordering along this line but there is a criticality in the 
paraphase. The susceptibility above the zero-temperature critical (T c = J) is 

t Q = 2e~ 2J/T . (70) 

This exponential behaviour corresponds to the critical exponents 7 = v = 00. 
Line aa: Along this line 

t - 4J|T T 7 Tcl e- 2r ^ (71) 

and, therefore, for both T <T C and T > T c , the susceptibility x ~ l/l^ol obeys the scaling 
law x — Xo/\T — T c \ with an exponentially increasing scaling amplitude xo- There is a 
smooth crossover between the pure exponential behaviour (70) and the power law (71). As 
(u /T) w J 2 /8r 3 for T ~ J > T, we have \(j) \ ~ \t \^ 2 , where t is given by Eq. (71). The 
order parameter exponent (3 has the classical value (3 — 1/2 but the scaling amplitude 
[~ T c _1 exp(— r/T c )] of the order parameter O exponentially decreases for T c — > 0; for 
a comparison of these low-temperature MF properties with high-temperature ones; see 
Ref. dHDj. 

The Eqs. (7), (70) and (71) yield a crossover exponent u = 00. Perhaps, this value of vq 
will be different provided the fluctuation phenomena are taken into account. 

6. 3. Lowest order perturbation theory 
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The (T, T) domains of validity of the MF results in Sec. 6.2 can be investigated by the 
Ginzburg criterion; see, e.g., Ref. pQ. Here we shall use the standard derivation of this 
criterion from the first-order perturbation contribution to the "self-energy" to: 



to = t + UuoA^to) 



(72) 



with 



Ai(t Q ) = 



d d k 
(2n) d 



Si(t ,k) , 



(73) 



where 



Si(t ,k) 



r cth[(r/T)VF+^ 



o 



(74) 



Neglecting a term of order Uq we can substitute to in the intergal Ai(t ) by t . The 
Ginzburg criterion for the validity of the MF results will be given in the general form 



The same criterion defines the validity of the Gaussian approximation of noninteracting 
fluctuations above and below the critical point. The general criterion (75) is valid for both 
the paramagnetic (t > 0) and the ferromagnetic (to < 0) phases; in the ferromagnetic 
phase a factor 1/2 should be added to the r.h.s. of (75). 

Obviously the results from this first-order perturbation approximation depend on the 
properties of the integral (73). We shall write this integral in the form 



where y = Xk, and g is given by Eq. (3). For high temperatures (A < 1, g < 1) the 
main contributions in the integral are from the small wave numbers y = Xk <C 1, and the 
approximation cthy ~ 1/y yields the standard perturbation integral A\ known from the 
classical theory (the Curie- Weiss limit). One can rederive the classical theory (cjj = 0) 
by using the high-temperature value of the interaction parameter: uq = (J 2 /12T 2 ); see 
Eq. (64) for T < T. 

In the low-temperature range of temperatures (A ^> 1, AA ~ 1) we can distinguish between 
the quantum limit g ^> 1 and the classical limit g <C 1. In the quantum limit (g 3> 1) 
the coth can be approximated with unity by neglecting exponentially small correction 
terms. They are very similar to those given by Eq. (69)-(71). These corrections enter in 
the renormalized parameter t and give small (~ w ) corrections to the coefficients of the 
corresponding exponential terms in t ; see Eqs. (70) and (71). 

Let us denote the integral which is obtained from A\ for cthy ~ 1 by Aoi. The difference 
uq[Ai — Aqi\ has been estimated [IHj to be of order A -2 ~ T 2 for 3<i-systems; note that 



to | > 12«o[A x (0) -A 1 (t )} . 



(75) 




(76) 
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in the low-temperature limit, u ~ T. This type of temperature corrections to the pure 
quantum limit has been widely used in interpretations of experimental results for quantum 
ferroelectrics [2T H I15fij . The same corrections were also derived in Refs. |154| 1155] : for a 
calculation in the Hartree limit, see Ref. 161J . 

At extremely low but finite temperature (T ~ 0) one may try to take into account the 
temperature corrections to the critical behaviour at the absolute zero (T = 0) as an 
alternative to the present consideration, where we take the low temperature limit in the 
general self-energy integral (76). Taking into account the temperature corrections to the 
zero-temperature limiting case is equivalent to an estimation of the difference arising 
from the substitution of summation over the Matsubara frequencies with an integration 
according to the rule (35). It is clear that such corrections will come from the counterterm 
series in powers of g~ 2 in the Euler-Maclaurin summation formula. Therefore, these 
corrections could not be given by d— dependent powers in T as is in the Hartree limit 
considered in Ref. |161| . 

The temperature corrections are important for the phase transition properties along lines 
like AJ and aa, where r ~ J. Within the present lowest order perturbation theory the 
corrections in powers of T should be considered small. However, it has been proven by 
the "parquet" summation ^H] that they are quite big and essentially influence the critical 
behaviour in the classical low-temperature region (g < 1, A ^> 1). 

Now we shall consider the Ginzburg criterion (75). In the quantum limit we must sub- 
stitute the integral A\ with A m . As a result of CQC, the integral A i yields the upper 
and lower borderline dimensionalities: dy = 3 and dy — 1. In order to simplify the cal- 
culations we shall consider the case d = 2 which shows the main features of the quantum 
critical behaviour for all dimensionalities 1 < d < 3. 

The straightforward calculation gives the criterion (75) in a simple form: 



This criterion cannot be applied to the T-transitions along the lines A J and aa where the 
behaviour is classical (A < £). The reason is that the strong quantum condition (g 3> 1) 
has been used in the derivation of inequality (77). 

Along the line bb we obtain the criterion 



If we set in (78) r c = J, we shall find a criterion along the line 0J. Because of the very 
sharp slope of the transition curve T C (T) near T = J, the quantum criterions along the 
lines bb and 0J are practically the same. The Ginzburg critical region (~ 10 _2 r c ) given 
by (78) is well established. It enlarges at T > by T-correction terms. 




(77) 




(78) 
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The quantum condition (3) along the line bb can be written in the simple form J(|T — 
r c |) 3> T 2 provided (T/r c ) 2 ps 1. As J ^> T, this condition becomes 

|r-r c |>T, (79) 

which is obviously consistent with (78). The condition (79) shows that the quantum 
region approaches the critical point T c when T decreases to zero. The whole surrounding 
of the zero-temperature critical point (T c = J) is influenced by quantum effects which 
produce the quantum fluctuation region. 

A more detailed analysis including the onset of thermal fluctuations can be obtained by 
the exact calculation of the integral A\ at d — 2: 



= 1 ln shVM+Z 

The approximation shy ~ exp(y)/2 corresponds to Ai ~ Aqi. The thermal corrections 
are obtained in powers of (g 2 /\A) which, together with the factor w ~ T, yields the 
lowest order correction term to Wo^oi °f the type ^ 2 |^o|- 

The asymptotic critical behaviour corresponding to the T— transitions will not exhibit 
quantum critical phenomena but rather a low-temperature classical behaviour which is 
different from the high-temperature one. The only quantum effect is related with CQC 
at zero-temperature T— transitions (Sec. 6.4). A study of the Ginzburg critical region 
for zero-temperature structural phase transitions has been presented in Ref. |162j . The 
results are consistent with the present analysis. 

6.4- Renormalization group arguments 

The RG recursion relation for the interaction parameter w in the one-loop approximation 
will be 

u' = b 4 - d u [l - 36u A 2 (0, b)} , (81) 

where 

Mi»i) = (82) 

In Eq. (82) Ai(to,b) is the integral (73) with lower (0 < b^ 1 < 1) and upper (A = 1) 
cutoffs of the wave number k. The other recursion relations are |ll()j 

t' = b 4 ~ d [t + UuoA^to, b)} , (83) 

and 

A' = b- z X , (84) 

where the dynamical critical exponent z is equal to unity in this order of the theory. 

In the high-temperature region (A < a) the integrals Ai and A 2 can be substituted with 
the classical integrals by setting cth(y) ~ 1/y. In this case the integral A 2 (0, b) has a 
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logarithmic infrared divergence at the upper critical dimensionality dy = 4. For d = 4, 
A 2 (0,b) =K d \nb. 

Further, using the standard RG analysis [20] one reveals the usual universality class of the 
critical behaviour of the classical Ising model (r = 0). Besides, there is a possibility to 
perform a calculation of dynamical critical exponent z and, hence, to reveal the quantum 
dynamics of TIM |163j . The e corrections to z are calculated from the g-dependent self- 
energy diagrams in two- and higher-loop approximations. These corrections are small 
compared with the e-corrections to the static exponents, for example, v and 7. 

The low temperature limit is treated as shown in Sec. 6.2. In this case we should make 
the approximation cthy ~ 1 in Eqs. (81)-(83). Thus we recover the dimensional CQC: 
d — > (d + 1) = D. The analysis in e = (3 — d) yields results for the quantum critical 
behaviour corresponding to the nontrivial Ising universality class (-0,1) (2D]- For all 
dimensionalities d > 3 the QC behaviour will be described by the classical MF universality 
class, while for dimensionalities 1 < d < 3 the QC behaviour will be nontrivial. This 
result is straightforwardly generalized for a n-component real field with a bare correlation 
function of the form (62). 

Thus we obtain that in TIM CQC satisfies the universality property discussed in Sec. 4.2. 
An important point in the mechanism of this crossover is that the interaction parameter 
u in the recursion relation (81) is changed to the parameter v = (u /X): 

V ' Q = b 3 ~ d v {l - 36v K 3 lnb] . (85) 

The quantum effects play the crucial role for CQC in TIM. If the RG equations are treated 
in the classical scheme, as shown in Refs. [1531 1155] . the transformation of Eq. (81) for uq 
to that for vq cannot be performed in the way shown by Eq. (85). An extra-factor A -1 
that remains in the second term of Eq. (85) leads to the prediction of a Gaussian critical 
behaviour at low temperatures for dimensionalities d > 2. 

The RG investigation |1541 1155j of the critical behaviour along the onset T— transition at 
T c = yields critical exponents 7 = 2u = 2 instead of classical exponents 7 = 2v = 1 
(here we ignore e— corrections). These exponents do not depend on the dimensionality d 
as is in the Hartree limit [161J. For to — 0, i.e., on the line AJ, the integral A\ has a 
logarithmic divergence which is an evidence of strong thermal fluctuations. Perhaps this 
nonuniversal critical exponetns of the zero temperature T— transition are the outcome of 
the simultaneous effect of classical and quantum fluctuations. 

The HLTC of the T-transitions from classical critical exponents at high temperatures to 
low-temperature critical exponents at the zero temperature phase transition (T c = 0), 
7 = 2i> = 2, shows that zuq = 1. Therefore, within this accuracy of the theory one cannot 
conclude whether the cirterion (9) is satisfied or not. These bare values of z and vq have e 
corrections from RG and we should have in mind that the correction to z is of order 0(e 2 ) 
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whereas that to v is of first order in e with a positive sign. Therefore the criterion (9) is 
not satisfied and, hence, the asymptotic critical behaviour along the onset T— transition 
to the absolute zero (T c = 0) is classical. 

Another important feature of the quantum critical behaviour is that it is unstable with 
respect to any perturbation of the temperature from zero. This feature is common to the 
quantum critical phenomena in all systems. Of course, the quantum critical phenomena 
will belong to the nontrivial ((i-dependent) class of critical behaviour in the quantum 
critical region shown by Eq. (78). The zero-temperature transitions outside this narrow 
region will exhibit the usual MF behaviour. It becomes clear from this picture that like in 
the XY model discussed in Sec. 5.6, the uniaxial ferroelectics and ferromagnets described 
by TIM, are convenient for the experimental observation of quantum critical phenomena 
near T-transitions at low and extremely low temperatures. 

7. Disorder effects 

Here we shall review investigations on the effects of quenched disorder of type "random 
critical temperature" which is caused by randomly distributed (quenched) impurities and 
inhomogeneities [H H21 I164[ I165[ I166j . The main results [1683 H69j for the quantum 
critical behaviour in systems with disorder of type "random field" will also be mentioned 
(for the random field disorder see, e.g., Refs. [Tj 11641 fl65j ) . 

The critical behaviour in disordered systems obeys the Harris criterion derived by rather 
general arguments |167j . This criterion states that the disorder will be irrelevant to the 
(classical) criticality, if the following inequality between the specific heat exponent a and 
the correlation length exponent v takes place: [av — 2) > 0. One can check that this 
criterion is satisfied in all examples considered below. 

The first RG investigation |17U| I171j of the quantum critical behaviour in disordered 
quantum systems reveals an instability of the usual quantum critical phenomena at T = 
towards the disorder of type random impurities. The instability has been proven |171j for 
a quite general quantum field model which describes the quantum phase transitions in 
almost all known quantum systems. This problem is closely related with the description of 
the localization in real superfluids. Here we shall consider the instability phenomenon, RG 
methods of stabilization of the quantum criticality and the related problems of localization 
and superfluid-insulator transition as well as the formation of "glass-like" phases (2H I148j 
in disordered Bose fluids. 

1.1. Random impurities 

The disorder of type random impurities changes the local interaction responsible for the 
phase transition and, hence, the critical temperature which depends on the sites of the 
crystal lattice. In the continuum limit the local (nonequilibrium |2()j ) critical temperature 
depends on the spatial vector r. This fact is taken into account in the effective Hamiltonian 
by an additional 2 -term containing a random function (p(r) that obeys the Gaussian 
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distribution 

g(f,f') = [(p a (r)ip a ,(r')] R = A6 aa >6(r-r') , (86) 

or, in the fc-space, 

g(k,k') = & a (k)p a/ (k')] R = A8 aa ,5(k + k',0) . (87) 

Here [ } R denotes the operation of averaging (hereafter, the suffix "R" will stand for 
quantities related to disordered systems). We shall often use the notations g aa i = 5 aa 'g 
and g aa > = 5 aal g. 

The distribution function (86) or respectively, g(r,r'), which is often called "random cor- 
relation function," describes quenched impurities with the so-called short-range random 
correlations. Sometimes, this type of random impurities is referred to as correlations" 
or, even, "5— impurities" because of their mathematical expression by the 5— function 
valid in the continuum limit. 

For the long-range random correlations, the 5-function in Eq. (86) should be substituted 
by a function of type f(R) = 1/R @ , where R = \R\ = \r — r'\, < 6 < d. In this case, 
the random distribution (correlation) function is given by g(r,r') = g(R) and 180J 

g(R) = ■ (88) 

A justification of this definition of long-range correlations is presented in an investiga- 
tion |168j of the random-field problem |T| 11641 [165 . The Fourier transformations of 
various shapes of random correlation functions g(R) with powerwise or more complex 
dependence on the distance R are considered in details, and the Fourier transforms 
g(k,k') = S(k + k')g(k) are given for various shapes of g(R) |168t 1169] . Note, that 
in the general case of random correlations, different from the 5— correlations given by 
Eq. (86), the parameters A and A in Eqs. (86) and (87) can be different from each other 
but both of them should be nonnegative on account of the requirement for the stability 
of the Gaussian distribution. 

An alternative way of representation |168| I169j of the correlation function g(R) from 
Eq. (88) seems to be more convenient for several considerations (see Sec. 7.4.1). This 
representation is given by the exponent 0' = (d — O), d > & > and is suitable for the 
scaling and RG analysis in the k— space. 

In order to describe the disorder of type random impurities a new term should be added 
to the quantum Hamiltonian (27) in the form 

n R [<t>\ = ^= p(*i-*2)&K*O0«K*2) . (89) 

This term can be derived by the Hubbard-Stratonovich transformation from the micro- 
scopic model of the concrete system or, alternatively, postulated and used in further 
investigations of disorder effects. 
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The sum of the Hamiltonian parts (27) and (89), i.e., the total Hamiltonian 

ft to t(0)= HW+HM)* ( 9 °) 
describes a great amount of quantum systems with quenched (randomly distributed) impu- 
rities, often called "impure systems." As usual, we shall discuss also the critical behaviour 
of the corresponding classical systems {oji = 0). The respective pure systems are those, 
for which the Hamiltonian part 7Yr(0) is equal to zero [ip(k) = 0, or, g(R) = 0, which 
gives the same results for the thermodynamic and correlation properties]. 

The thermodynamic behaviour strongly depends on the properties of the random po- 
tential represented by the random function ip(r). The random potential <p(f) is not a 
thermodynamic variable, that is why, the theoretical treatment will include two steps. At 
first one should try to calculate the thermal averages as functionals of (f(r). As a second 
step the averaging over the random function is done with the help of Eqs. (86)-(88). As 
a result one obtains the thermodynamic quantities as functions of the disorder parameter 
A; see. e.g., pTI IT21 \TM ITES] . 

Alternatively, an ansatz called the "replica-trick" can be used (see, e.g., |164| 11651 1166] ). 
In perturbative and perturbative-like investigations as, for example RG, this trick should 
always give the same results as the method of direct perturbation treatment of disorder. 
This fact has a simple explanation. The replica index (m) appears in the perturbation 
terms through polynomials and, hence, there is no problem in the analytical continuation 
of the results in the limit m — > 0. The limit must be taken at the end of the replica-trick 
calculation in order to retrieve the original disordered system |164[I165]I166| . In some non- 
perturbative calculations, however, the same analytical continuation and the justification 
of the limit m — > are not very easy; see, e.g., Ref. |164( I165[ 1166] . A recent theory of 
of the thermodynamic behaviour of disordered superfluids below the finite-temperature 
critical point (T c > 0) has been developed by Lopatin and Vinokur |172j on the basis of 
the Beliaev theory |173| I174j and the standard replica trick treatment of 5— correlated 
random impurities. 

The results reviewed in the remainder of this paper are obtained by RG and provide an 
information about the asymptotic critical behaviour, namely, the critical behaviour in an 
infinitesimally close vicinity of the critical point. Note, that less information is available 
about the MF behaviour of systems with quenched disorder and the pre-asymptotic critical 
behaviour. 

7.2. Results for classical and quantum models at finite temperatures 

The RG investigation of quantum effective Hamiltonians (27)-(29) with the additional 
term (88) for all possible values of m, ml and a was published for the first time in Refs. |17U] 
I171j . Here we shall review the results from Refs. |170] I17T] together with the necessary 
information from other papers; see, e.g., Refs. [T2" | 1164] 1165] 1166 . 

7.2.1. Some results for classical systems 
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In systems with symmetry index n > 4 the disorder effects are irrelevant to the critical 
behaviour at finite temperature critical points (T c > 0). The finite temperature critical 
behaviour of systems with symmetry indices n > 4 is described by the usual Heisenberg 
FP of the corresponding pure systems (A = 0) (for n = 1 this FP is usually called Ising 
FP) [HHSUHj. Here we shall often refer to this FP as "pure" FP or, shortly, PFP and 
even P. Note, that PFP exists also for some values n < 4 but for such systems this FP is 
unstable towards disorder (A > 0). 

For systems with n < 4 the disorder effects essentially influence the critical behaviour at 
finite critical points (T c > 0); see. e.g., Ref. [T2]. In this case the critical behaviour is 
described by a "random" FP established for the first time by Lubensky |175j : hereafter 
referred to as RFP, or Lubensky RFP, or shortly denoted by "R" . 

The effect of the random impurities on the critical behaviour in Ising-like systems (n = 1) 
is quite peculiar (see, e.g., Refs. [T2J I165j ). The RG equations exhibit a special form 
of degeneration and, as a consequence, the usual perturbation expansion breaks down 
and should be conveniently modified. For example, the usual e— expansion should be 
substituted with an expansion in non-integer powers (~ e P//2 , p — 1, 2, ...) of e = (du — d). 

The modified e— expansion yields special FP which attracts the RG flows through spiral 
paths. Such FPs are often called "focuses" " or focal FPs. These FPs give complex values 
for the correction-to-scaling exponents J2| which result in an oscillatory behaviour of the 
corrections to the scaling laws. Perhaps, this extraordinary aspect of the impure critical 
behaviour is a feature of systems where an anisotropic ordering appears below the critical 
temperature or other anisotropy effects are present, as shown in Ref. [HZ| for a wide class 
of anisotropic impure systems. We shall not dwell on the case n = 1 (for this case, see, 
e.g., Refs. ^]) rather we shall discuss systems with continuous symmetry (n > 1). 
Let us, however, mention that several FPs of focal type are discussed in Sees. 7.4 and 
7.5 for other reasons which are associated with the effective increase of the upper critical 
dimensionality d§ of the disordered system (d§ > du)- 

7. 2. 2. Quantum systems at finite temperature 

The picture outlined in Sec. 7.2.1 for the asymptotic critical behaviour of classical systems 
does not change for quantum systems at finite critical points (T c > 0). The additional 
information which can be obtained for the finite temperature critical behaviour of disor- 
dered quantum systems is about the quantum critical dynamics, in particular, about the 
value of dynamical critical exponent z. The results |170| I171j for the dynamical critical 
exponent z of finite temperature critical points can be written in the form: 

(A — n) 

Z = a+ T< TV 6 ' 1<™<4, (91) 

8(n — 1) 

where e = (4 — d). Eq. (91) has been obtained for the XY model and Bose systems 
which are described by the bare correlation function (28). In the case of models that have 
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Hamiltonians (27) with the correlation function (29), the dynamical exponent has been 
obtained |17()[ I171j in the form 



z = — + 



a (4 - n) 



m=0 



1 < n < 4 . 



(92) 



m 8(n — 1) 



For m' > the exponent z is given by Eq. (36). In the latter case the critical exponent z 
has no correction to the first order in e |171j . These results correspond to RFP, namely, 
when the disorder is relevant. The same results |171j have been confirmed for disordered 
itinerant antiferromagnets [17(11 H77| 

Eq. (91) coincides with that obtained in Ref. |178j for disordered classical systems within 
an approach based on a time dependent Landau-Ginzburg equation [T2J 122] • The results 
(91) and (92) demonstrate that the dynamical critical exponent z in systems with random 
impurities has an e— correction in the one-loop (first order in e = 4 — d) approximation 
whereas the dynamical exponent of the corresponding pure system (A = 0) has an e- 
correction of order 0(e 2 ). This is a direct consequence of the fact that the random 
function (p(k) does not depend on the Matsubara frequency uoi. 

7. 3. Instability of the quantum critical behaviour in disordered systems 

The zero-temperature critical behaviour exhibits an instability with respect to the quenched 
disorder |17U| I17lj . The mechanism of the instability is related to CQC. Let us consider 
the RG equations corresponding to this problem |171j . In the limit T c — > the parameters 
of the Hamiltonian obey the lowest order RG transformations of the type |17Ut 1171] : 



where zq denotes the bare value z(0) of the dynamical exponent z(e). The "tree" approx- 
imation used in Eqs. (93) and (94) is sufficient for our present discussion (the one-loop 
RG equations have been derived in Ref. |171j ). 

Because of CQC the interaction parameter v is relevant to the zero-temperature critical 
behaviour in spatial dimensionalities d < (2a — zq), whereas the disorder parameter A 
is relevant to d < 2a. This is readily seen from Eqs. (93) and (94). The difference in 
the upper critical dimensionalities of the parametes v and A is produced by CQC and 
the lack of c^-dependence of the random function <p(k). On one side, at the classical 
borderline dimensionality djj = 2a, the parameter v is irrelevant and, hence, the RG 
equations in d = (4 — e) dimensions describe a simple Gaussian instability towards the 
disorder parameter A. On the other side, the zero-temperature e-expansion in terms of 
eo = {2a — d — zq) yields pure FPs of the type A* = that are unstable with respect to 
the disorder parameter A |171j . In these two variants of the theory FPs of corresponding 



v 




(93) 



and 



A' = b 2a - d A , 



(94) 
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Figure 5: A scheme representing four relevant FPs (P, R, Gl, Un) in disordered systems. 



pure system are unstable with respect to disorder effects for dimensionalities less than the 
upper borderline dimensionality du |171j . 

Moreover, the disorder itself does not generate new stable FPs in the quantum limit 
(T c — > 0). This is a clear indication for the lack of a standard (pure or random) zero- 
temperature critical behaviour and, therefore, one should expect either some unconven- 
tional (multi) critical behaviour in these systems at zero temperature or a first-order phase 
transition. 

The zero-temperature and, perhaps, low-temperature instability of the standard critical 
behaviour in pure systems against disorder effects of type quenched impurities is probably 
related with the instability of the usual ground state of the respective system and, hence, 
with localization effects in real disordered superfluids (see Sec. 7.7.3). Therefore, it seems 
important to have a good notion about the relevant FPs which govern the impure quantum 
systems. Using the results discussed up to now we can distinguish between four types of 
critical behaviour of pure and disordered Bose fluids and XY systems described by the 
Eqs. (27), (28), and (90). They are given by four FPs: (i) Pure FP (P) corresponding to 
T c > and A = (lack of disorder), (ii) Gaussian-like FP (Gl), corresponding to T c = 
and A = 0, Random FP (R) - for T c > and A > 0, and Unstable FP (Un) corresponding 
to the case of zero-temperature critical point (T c = 0) of disordered system A > 0. These 
main four FPs are shown in Fig. 5 at their places in the parameter space (A,T). 

For other quantum systems, for which Eq. (29) should be used instead of Eq. (28), the 
diagram in Fig. 5 remains valid, provided the Gl FP is substituted with the usual zero- 
temperature FP corresponding to CQC in pure systems (2D]- Thus we have two types of 
instability of the quantum critical behaviour towards random impurities: (a) instability of 
the quantum critical non-universal behaviour represented by G1FP, and (b) instability of 
the quantum critical universal behaviour described by the universality PFP (see Sec. 5.4). 

The effect of quenched (random) inpurities on the quantum tricritical behaviour has been 
also investigated |179j . The RG results demonstrate an instability of the pure quantum 



57 



critical behaviour towards the disorder and lack of any new stable FP of the RG equa- 
tions |179j . Therefore, the usual critical and tricritical points exhibit an instability due to 
CQC which takes place in the pure subsystem (the subsystem described by A = 0) but 
does not affect the "disorder subsystem," described by the Hamiltonian part (90). 

The random field disorder is itself a source of a dimensional crossover 11641 1165j . 

This is another type of quenched disorder that may occur in real experiments. It has been 
shown \WH\ I169| that the outcome of the simultaneous action of quantum and random 
field effects at zero temperature is again an instability of the quantum critical behaviour 
below the upper borderline dimensionality (the latter dimensionality is higher than 2a for 
the presence of short- or long-range correlations of the random fields |168| [L69 ) . 

7.4- Long-range random correlations and extended impurities 

7.4-1- Long-range random correlations. 

The quantum critical behaviour is unstable also towards disorder effects produced by 
random impurities with long-range random correlations considered for the first time by 
Weinrib and Halperin |18()j (for the calculation of the critical exponent r] in this case, see 
Ref. |181j ). The quenched disorder is described again by the Hamiltonian part (90) but 
the correlation function g(R) is given by Eq. (88). 

The Fourier transform of the random correlation function (88), to the leading order in 
the wave vector k, has the form g{k) = Ax > for > d, and 

g(k) = A 2 k @ ~ d (95) 

for G < d. These two cases have been combined in Ref. [1K01 . where the RG analysis is 
carried out by the correlation function 

g(k) = Ax + A 2 k @ ~ d , (96) 

which contains both short-range (A 2 = 0) and long-range (Ax = 0, G < d) random 
correlations between the quenched impurities. 

For classical models (ut = 0) [180J the RG analysis gives new stable FP of a focal type 
with complex eigenvalues, leading to oscillating corrections to scaling laws and consistent 
with the Harris criterion |167j . Obviously, the same features will remain valid for quantum 
systems at finite critical temperatures T c > 0. As there is no spatial anisotropy in this 
model and, moreover, no anisotropy is present at all in the n— dimensional space of the 
order parameter vector, we are faced with an example when focal FP, describing oscillating 
corrections to the scaling laws can appear also in cases of a total lack of anisotropy. The 
only reason, for which such a focal FP appears in this class of systems is the disorder-driven 
dimensional crossover associated with the upper borderline dimensionality (d§ > djj) as 
explained in more details at a next stage of our consideration. 
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Let us consider why the classical impure behaviour is proven to be stable by RG investi- 
gations, at least for a rather large class of systems |18()j . This is a result of an ansatz, by 
which a double e— expansion is introduced in the RG studies ;180j. However, without this 
ansatz, namely, by following the standard RG scheme, one should obtain an instability of 
the critical behaviour of the impure system as a result of a lack of any stable FPs of the 
RG equations in the nontrivial dimensionality range d^ < d < d§. 

The reason for the instability is in the dominant role of the disorder parameter A2 as- 
sociated with the long-range random correlations. The presence of long-range random 
correlations simply quells the the fluctuation interaction effects represented by the pa- 
rameter v as well as the short-range random correlations, represented by the parameter 
Ai. The behaviour of the system very much resembles the behaviour of noninteracting 
fluctuations (free field; v — 0) in a random potential. 

In order to clarify this problem one may perform a simple dimensional analysis jU of 
the Hamiltonian, or, equivalently, to perform the RG rescaling within the simple "tree" 
approximation. The result of the RG rescaling procedure is 

v ' = b 2a ~ d v, A x = b 2a ~ d ^ 1 , (97) 

and 

A 2 = 6 2ct -°A 2 . (98) 

If the difference (d — 0) > is of order unity, as is for real systems, an instability of the 
classical critical behaviour corresponding to the pure system (Ax = A 2 = 0) will occur 
and, moreover, a new stable critical behaviour will not appear at all. The reason is that 
the parameter A 2 will always be relevant and the perturbation contributions which are 
supposed to be of order e cannot compensate the growth of this parameter to infinity. 

The ansatz introduced for the first time in Ref. j!8Uj for this type of systems consists 
in the assumption that, at least during the formal derivation and analysis of the RG 
equations, (2a — 0) = 5 will be a small quantity of order e = (2a — d); remember, that 
the upper critical dimensionality of the pure system is du = 2a. Then one performs a 
double (e, S) which yields a new stable FP of focal type |18()j and complex exponents 
for the correction-to-scaling laws which describe the critical behaviour of systems with 
long-range random correlations. 

The formal assumption that 5 = (2a - 0) <C 1 makes possible to develop a double 
e— expansion but the form of Eq. (98) does not allow an easy way to obtain the upper 
borderline dimensionality d§ of the disordered system, which seems to be higher than 
that (du = 2a) of the pure system on account of the particular effect of the parameter 
A 2 . This is a peculiar feature of the above described approach because usually the tree 
approximation immediately yields the upper critical dimensionality. The reason is in 
the lack of a ci-dependence in the scaling transformation (98). Thus one is left with the 
opportunity to determine the dimensionality d§ from the singularities of the one-loop 
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integrals rather than from the lowest order tree approximation. Therefore, one begins 
with the derivation of the RG equations for a general dimensionality d and only after 
that one has to determine the critical dimensionality d§, at which RG equations are to 
be expanded. 

The formal comparison of the scaling (b—) factors in Eqs. (97) and (98) shows that 
(2a — 0) = [e — (d — 0)] and, hence, one may speculate that the upper dimensionality 
d§ of the disordered system is raised with respect to djj by (d — 0) but this argument, 
although correct (see below), cannot be used to determine d§ as a number independent 
of d. The reason for this difficulty is in the choice of the exponent that decribes the 
R— dependence in Eq. (88). If we use the exponent 0' = (d — 0), which has been already 
introduced in our analysis, the Eqs. (88), (95), (96), and (98) will change in an obvious 
way. For example, the k— dependence in the second term in the r.h.s. of Eq. (96) will be 
written as (l/k e> ), and the b— factor in Eq. (98) will take the form fr( 2<7 + e Now we 
can be certain that d§ is exactly equal to (2a + 0'). The difference (d\j — du) is equal to 
0' = (d— 0) = (5 — e) > 0. Note, that the focal FP revealed in Ref. |18Uj has a singularity 
just at 5 = e. 

The main result of the RG analysis within the double e— expansion - the stabilization of 
the quantum critical behaviour in disordered systems of special type, has been achieved 
by an ansatz which does not affect the (1/T)— axis of the total (r, x) space of the critical 
events. It is easy to see that CQC in the limit (T — > 0) will cause an instability of all 
available finite temperature FPs in these systems. 

7.4-2. Extended impurities 

Up to now we have discussed random "point" impurities, i.e. impurities localized at single 
points r. The disorder of type extended (non-point) random impurities is described by 
infinitely-ranged random correlations along one or more spatial dimensions, say, d < d, 
and short-range random correlations along the other (d — d) dimensions |18H| I184[ 11851 

HEEI- 

For example, in systems with a slab geometry the appropriate choice of extended impuri- 
ties is the one-dimensional (linear) impurities oriented along the direction of the small size 
L (the slab thickness) and randomly distributed along the other spatial directions |182j . 
The extended quenched impurities are described by a modification of the model for short- 
and long-range correlated point impurities, in which these finite-range correlations along 
the "small size" Lq, for instance, along the z— axis, are substituted with infinite-range 
correlations. The random correlations along the long dimensions of the slab, i.e., along the 
x— and y— axes, are kept of the short-range type. The length scale of "the infinite-range" 
correlations is much larger than the correlation length £ and the thickness Lq of the slab. 
So, the strongly correlated along the small size point impurities behave like continuous 
uniform strings. In regard to the critical behaviour this disorder acts like point impurities 
with a short-range random distribution along the large (infinite) dimensions L, and an 
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uniform distribution along the small size L [182 . 

In infinite systems one may consider: (i) straight lines of impurities or straight dislocation 
lines of random orientation, (ii) two-dimensional (plane, d = 2) defects, and finally, (iii) 
randomly distributed impurities or defects of dimensionality d < d. 

Now the important question is about the way of description of infinite-range correlations. 
This is simple because one takes the exponent of the correlation function (88) to tend to 
zero (G — > 0) along the directions in the d— dimensional subspace of extended impurities. 
Then the correlations along this directions will be of infinite range but uniform. The body 
will be homogeneous (aligned) along this directions and the disorder will be present only 
along the rest part of spatial directions (d — d). Now we can suppose either short or long 
range random correlations along these directions and let us choose, as in many papers, 
the former opportunity. Then we shall have no defects (ideal order) along chosen axes 
of number d and a short-range disorder of type random impurities along the rest (d — d) 
axes of the d— dimensional coordinate space. The respective correlation function of type 
(86) is given in the form [H3 1113 CH3 QB5| 

g(f,r') = A5 d - J (f-f') . (99) 

The case when d = 2 and d — 1 corresponds to the exactly solvable McCoy- Wu lattice 
model |187[ll88(ll89j of disorder along one of the spatial directions for the two-dimensional 
Ising model. We must emphasize that the mentioned correspondence is only in some 
general features. 

In all cases the extended impurities will generate a critical behaviour instability for di- 
mensionalities less than d < d§ = (2a + d). The mechanism of this instability lies in 
the appearance of the finite shift d = (d§ — du) of the upper borderline dimensionality, 
i.e., the reason is similar to that producing instability by long-range random correlations 
(Sec. 7.4.2). In fact, the raised dimensionality d§ corresponds to the disorder parameter 
A and the critical dimensionality associated with the fluctuation interaction v remains 
the same as in the pure system (du = 2a). This means that the disorder effect represented 
by the parameter A suppresses the fluctuation interaction ~ v <f) A in Eq. (27) just like in 
the case of long-range random correlations described in Sec. 7.4.2 and, as a result, stable 
FP is lacking and the critical behaviour is unstable. 

The instability is overcomed again by an ansatz of a suitably chosen double e— expansion [183 , 
I184[ H851 H86j , namely under the supposition that d is a small quantity d ~ e and can be 
treated on the same footing as the expansion parameter e = (2a — d). Thus using d = e<i as 
a second small parameter in the double (e, e^)— expansion, the RG theory |1831 1184[ I185j 
gives a stable finite temperature critical behaviour by stable FP of a focal type. This 
(Dorogovtsev) FP is different from focal FP discussed in Sec. 7.4.1 and describes the crit- 
ical behaviour with complex corrections corresponding to the extended impurities with 
a random distribution (99). Despite of the number of similarities between the critical 
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phenomena, from one hand, in systems with point impurities with long-range random 
correlations and, from the other hand, in systems with extended inpurities with short- 
range random correlations there are also essential differences. These differences are mainly 
in theoretical aspects, namely, in particular features which create methodical problems. 
From an experimental point of view, however, the effects of these types of disorder on the 
critical behaviour can hardly be distinguished. 

The stability of the critical behaviour in disordered systems with extended impurities 
vanishes when the quantum limit (T — > 0) is taken and CQC occurs. The reason is again 
in the appearance of CQC which produces an effective suppression of the fluctuation 
interaction represented by the parameter v. 

7.4-3. Summary 

We can conclude that the quantum critical behaviour is unstable against all relevant 
forms of disorder used in the theoretical physics for the description of real systems. The 
instability of the quantum critical behaviour in pure systems with respect to disorder 
effects of several types considered in Sees. 7.3 and 7.4 is a result of dimensional CQC 
that suppresses the thermal fluctuation interactions but does not affect the disorder. 
The reason is simply in the fact that the quenched disorder is considered as time (r— ) 
independent. 

The similarity between dimensional CQC and the dimensional shift (d§ — djj) of the 
upper borderline dimensionality established in certain types of disordered classical systems 
can be used in the development of a unified approach to the description of dimensional 
crossover phenomena in quantum systems. It has been pointed out |185j that the related 
with the time (r— ) extra dimensionality of the quantum Ising model can be supposed 
throughout a double (e, e^)— expansion equal to in a complete analogy with the case of 
linear random impurities in classical models. 

7.5. The stability ansatz about the quantum criticality in impure systems 

The idea of Boyanovsky and Cardy [185J mentioned at the end of Sec. 7.5, has been 
applied by Weichman and Kim |190j for solving the problem for the instability of the 
critical behaviour at T = 0. In Ref. |19Uj disordered superfluids with random impurities 
have been investigated. These systems are described by the sum the Hamiltonian parts 
(27) and (88), and the bare Green function, given by Eq. (28). But the results in Ref. 
can be applied beyond the model of disordered Bose fluids and include other quantum 
systems. Note, that an artificial upper cutoff A w for the Matsubata frequencies uii has 
been used in this investigation |19Uj but this cutoff is avoided in a later paper |115j . 

Weichman and Kim |19()j considered linear impurities along the time (r— ) axis instead 
of the spatial ones. Thus the impure strings along the time direction are of size approx- 
imately given by f3 — 1/T. For A ~ 1/vf — > (for simplicity, here we set 9 = 1/2), 
when the temperature time (r— ) dimension reduces to zero and one can take the contin- 
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uum limit along the axis", these lines shrink to randomly distributed points in the 
d-dimensional volume of the system. For T — > 0, however, the "length" f3 of these lines 
tends to infinity. Note, that the consideration based on the exponent 8 = 1/2 corresponds 
to a dynamical exponent z = 2 and, hence, to all systems described by the Hamiltonian 
(27) with exponents a, m', and m, that satisfy Eq. (36) for z = 2. But the genaral features 
of this considerations are valid for all systems described by Eqs. (27), (28), and (29) and 
the respective disordered systems given by the Hamiltonian (90). 

In fact, within the framework of the problem, considered in Ref. |190j and later in 
Ref. |115j . the dynamical critical exponent z is equal to two and, hence, the term "lines" 
should be used, if we consider a (d + l)-dimensional space-time(r). However, if we con- 
sider this space-time representation equivalent to (d + 2) dimensional space, as is in CQC 
dimensional crossover for Bose fluids, we should stretch the "lines" (or strings) of "length" 
(3 ~ A 2 in planes of area A x A. In general, we shall have hyperp lanes of dimensionality 
z and "area" ~ X z . For A — > the randomly distributed objects will shrink to randomly 
distributed point impurities. These are two different geometrical representations of the 
same phenomenon. The formal analysis can be performed with the additional small pa- 
rameter e' d which runs from zero to unity and thus describes the classical (e' d = 0) and 
the quantum (e' d = z) limits of Bose system with point impurities. Alternatively one may 
consider the small parameter = e'J z |190[ 1115] . The double e— expansion of the RG 
equations again yields stable focal FP. 

In result, the following interesting picture of the Bose systems critical behaviour is ob- 
tained |19()j . In the high-temperature range, where X < a, the critical behaviour is 
described by the high-temperature RFP [175 J corresponding to short-range (point-like) 
impurities, i.e, Lubensky RFP (see also our notations defined in Sec. 7.2.1). In the zero- 
temperature limit (T — > 0), the impure critical behaviour described by Lubensky RFP un- 
dergoes a crossover to a critical behaviour governed by extended (plane) impurities. This 
zero-temperature critical behaviour is described by the corresponding zero-temperature 
variant of Dorogovtsev RFP with complex stability exponents |183[I184| . In the Hamilto- 
nian parameter space (fi,v,A), Lubensky RFP is conjugate to Heisenberg PFP (A* = 0) 
and both of them are far from the zero-temperature plane f3 = oo. Dorogovtsev RFP is 
conjugate to G1FP (A Gi = T Gl = 0) considered in details in Section 5.3. Like Heisenberg 
PFP, G1FP is stable only for pure Bose fluids (for real bosons, n = 2). For any A > 
and T > 0, the RG flows tend to the finite-temperature Lubensky RFP, whereas the 
zero-temperature Dorogovtsev RFP will be attainable, if only the system is at T = 0. 

These four conjugate FPs describe the critical phenomena in disordered superfluids and 
all disordered quantum systems presented by the Hamiltonian given by Eqs. (27), (28), 
and (90). One may visualize this scheme by a picture similar to that given in Fig. 5. 

For the rest part of disordered quantum systems, namely, those described by the bare 
Green function from Eq. (29), the same scheme works again provided one substitutes 
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the nonniversality G1FP with the respective universality PFP (see also Sec. 5.4). Then 
the zero-temperature pure critical behaviour obeys the exact universality rule (Sec. 4), 
that is, the <i-dimensional zero-temperature critical behaviour is described exactly by 
the finite-temperature universality class (d + z,n), to which the system belongs (see also 
Sec. 5.4). The same is valid for the respective disordered systems. In this case, the 
disordered quantum system enters in the high-temperature universality class (d + z, n) of 
the respective classical system (ui = 0). 

7.6. Related problems 

7.6.1. Disordered ferro- and antiferromagnets 

An approach identical to the Weichman-Kim approach |19Uj (Sec. 7.5) has been used in 
an investigation of the quantum phase transition in disordered itinerant quantum antifer- 
romagnets \17ti\ 1177] . In the notations given by Eqs. (27)-(28) these systems correspond 
to Eq. (29) with m' — 0, m— 1. Note, that in Ref. |176j the result (92) for the dynamical 
critical exponent z has been generalized with the help of (e, ed)-expansion. 

The picture of the quantum phase transition properties of pure and disordered itinerant 
ferromagnets given in the theoretical treatment |191[ 11921 H93j is worth mentioning. The 
energy spectrum e(k) of the fluctuation mode (the magnetization) in the effective field 
Hamiltonian depends on the spatial dimensionality d: e(k) ~ k^ d ~ 2 \ The quenched 
disorder of type random impurities produces a diffusive electron dynamics which induces 
an effective long-range electron spin interaction of the form 1/i? 2 ^ 1- ^ and this gives a 
d— dependent spectrum of the type l/k( 2 ~ d \ The bare Green function G(q) of the same 
magnetization fluctuating modes is described with the help of the terms given by Eq. (29) 
for m! = 2, m — 1, and a = min(2, d—2) and an additional term of type k 2 which is often 
redundant. The d— dependence of the energy spectrum of the fluctuating order parameter 
leads to a completely new, and in some aspect, quite non- universal picture of the critical 
behaviour in the most interesting domain of spatial dimensionalities (d < A). 

In order to evaluate this result one should be acquainted with the behaviour of the respec- 
tive pure system, where the disorder is not present. By the same method of treatment - a 
derivation of an effective field Hamiltonian from the microscopic electronic Hamiltonian 
with the help of the Hubbard- Stratonovich transformation and by a generalization of the 
Hertz [20] approach the respective energy spectrum of the critical magnetic modes in pure 
itinerant ferromagnets was found again to be d— dependent [193] . Now the bare Green 
function Go(q) will be given by Eq. (29), provided the exponents have the values m! = 1, 
m = 1, a = 2 and, moreover, an additional term of c'k^ d ~^ with d > is added to the 
r.h.s. of the same expression for Go(q). As a result, a new universality class of critical 
behaviour arises at dimensionalities d < 3. This problem has a further development in 
Ref. |194j with the conclusion for a fluctuation-induced first order phase transition at 
dimensionalities d < 3. 
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For details about the motivation and the results of this challenging direction of research, 
see the original papers |191| I192| I193j . One should have in mind that the Hubbard- 
Stratonovich transformations as a general method of establishing a correspondence be- 
tween a microscopically formulated Hamiltonian and its effective field (quasi-macroscopic) 
counterpart have shortcomings unless one is interested only on the asymptotical long- 
wavelength limit (see, e.g., Refs. ["""""[ I195j ). On the other hand let us mention that the 
integration out of critical or auxiliary-to-critical modes in complex models almost always 
generates d— dependent terms in the spectrum of the principal fluctuating mode describ- 
ing the phase transition of interest to the particular study. In usual cases of effective field 
Hamiltonians with more than one fluctuating fields this procedure of a further reduction 
of the description is rarely correct. 

A density matrix RG theory |196] was applied to the one-dimensional (chain) random- 
exchange spin- 1/2 XXZ model |197j . The interplay of quantum fluctuations and disorder 
results in a disorder-induced quantum phase transition that exhibits a nonuniversal be- 
haviour of the spin correlations. 

7.6.2. Disorder in quantum Hall liquids 

A RG investigation of the effect of point random impurities with short-range correlations 
of FQHE systems was done by Schakel j!15j (see, also, Ref. |lj). The effective field CSGL 
theory, discussed in Sec. 5.7.2, has been considered in Ref. |115j in a reduced variant, in 
which the gauge fields are omitted and the Bose field <f)(q), which interacts only with the 
random mode <p(k) subjected to 5— correlations is described by Eqs. (27), (28), (86), (89) 
and (90). In this case one retrieves the instability |171j (Sec. 7.3) and in order to ensure a 
derealization transition, one is forced to look for a stabilization of the ground state. Note, 
that the gauge fields (the Chern-Simons field and the vector potential of the magnetic 
field) cause other fluctuation effects (see, e.g., Ref. [Tl [2021 12*03*] ) that additionally lead to 
an instability of RG FPs. Thus one cannot easily distinguish between the effects producing 
the instability of the quantum critical behaviour within the framework of the complete 
CSGL theory. Note, that the reduction to a 4 theory makes possible the investigation 
of the net effect of the disorder on FQHE. 

In order to restore the stability of the quantum critical behaviour in this disordered 
system, Schakel |115j has proposed a method which departs from the Weichman-Kim ap- 
proach |190j (see also Sec. 7.5). The main disadvantage of the Weichman-Kim treatment, 
as mentioned in Ref. |115j . is the introduced by these authors an upper cutoff for the 
Matsubara frequencies. Note, that such a cutoff, after the paper of Hertz [2*0] . is used 
by a number of authors. On the one hand, this cutoff is considered as irrelevant for the 
asymptotic critical behaviour which is governed by low-frequency fluctuation modes, but 
on the other hand, as mentioned in Ref. |115j . the same cutoff is difficult to justify as 
it would imply a notion for a "discrete" rather than continuous imaginary time (r— ) 
variable. 
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In the Schakel approach the frequency cutoff is avoided and the convergence of the fre- 
quency sums (35) in the perturbation terms is saved by introducing a noninteger imaginary 
time (t— ) dimensionality e^: < e& < 1. The value ea = 1 describes the quantum limit 
(T — » 0), when the imaginary time r G [0, 1/T] runs over the values in the whole in- 
terval [0, oo]. In this limit one can apply the standard rule (35). When = we deal 
with the respective classical system where the temperature is sufficiently high and the 
r— dependence of the fluctuation field <f>(x) can be safely ignored. The intermediate val- 
ues < €d < 1 simulate a form of an approximate interpolation between the two limiting 
cases. 

The respective double (e, e^)— expansion reveals stable FP corresponding to the quantum 
regime (T = 0) of the disordered system (for details, see Ref. |115p . In some aspects this 
FP is similar to that found by Weichman and Kim |19()| . 

The Schakel approach works also for pure systems. In this case, the variations of the di- 
mensionality €d can be used to describe the classical-to-quantum crossover of the critical 
behaviour. Note, that this interpolation procedure is rather approximate, and the de- 
scription is exact only for the asymptotic cases (e^ = 0, 1). The approximation is a result 
of the substitution of the frequency sums in the perturbative terms by e^— dimensional 
frequency integrations. This problem is important for all studies based on double e— ex- 
pansions and has been investigated in details in Ref. j!82j (see also Sec. 7.6.3). Perhaps, 
this approximation is one of the reasons for the singularity separating the quantum and 
classical regimes within the Schakel approach |115j . 

7. 6. 3. Thin films with quenched impurities 

The formal conformity between CQC in quantum systems and the finite size crossover 
in systems of slab geometry |198l 1199} I2UU| !2Ulj mentioned in Sec. 4.3 as well as ideas 
discussed in Sees. 7.4-7.5 have been used in Ref. |182j in the study of disordered thin films 
with quenched point and extended (line) impurities. The method of work is similar to the 
Weichman-Kim [190 3 and Schakel |115j methods. The investigation has been carried out 
by a double e— expansion and reveals stable FPs which describe the critical behaviour of 
the disordered films. 

The second small parameter in the double e— expansion is the quantity 5 G [0, 1] which 
comes from an ansatz, discussed also in Sec. 7.6.2, and is used to define a substitution 
of lattice sums over the wave vector components with a 5— fold integrals. The ansatz is 
not new and defines a 5— integration that has two slighty different variants as shown in 

Ref. jug. 

The method of integration over wave vectors with noninteger dimensionality lies in the 
basis of all double e— expansions discussed in this review. In Ref. |182j the two small 
parameters e and S have been chosen in the form: e = (4 — D e ^) and 5 = (-D c ff — d) , where 
djj = 4, D e f[ = (d + S) is the effective (not the physical) spatial dimensionality of the slab 
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(film). Obviously, the quantuty D cS should depend on the ratio y = (L /0> where L 
is the thickness of the film. This dependence can be ascribed to 5: S(y). Thus the RG 
treatment in Ref. |182j has been performed by a double (e, 5)— expansion, where the small 
parameters are defined as deviations from the effective dimensionality D e g of the system. 

As shown in Ref. |182j the 5— integration as a method of calculation of perturbation terms 
introduces a systematic error of a finite magnitude, except for the limiting cases 5 = 0, 1. 
While the RG studies are related with infrared divergences this finite error produces 
only incorrect values of the FP coordinates but does not affect the physically measurable 
quantities as the critical exponents. Since the precise location of the FP coordinates is 
important for specific intermediate RG calculations but is not important for the physical 
predictions, the RG results obtained so far by double e— expansions are reliable. 

Furthermore, one may try to find a functional relation between the noninteger dimen- 
sionality 5— and the ratio y = (L /£). Obviously, the dependence 5(y) should have the 
properties: 5(y) —>■ for y — > 0, and 5(y) — > 1 for y — > oo. The method used in Ref. }182j 
confirms this properties but the attempt to obtain an analytical dependence 8{y) have 
not been successful. Moreover, the method does not provide a sufficient accuracy of the 
results for the intermediate values y ~ 1 correspondong to 5 ~ 0.5. 

The comparison of the RG results for disordered thin films |182j with the respective 
results for models of disordered quantum systems shows that there should be a formal 
equivalence between the universality features of the critical behaviour in thin films and 
that in quantum systems such as ferroelectrics and magnets described by TIM; see Eq. (29) 
for a = 2, m = 2, and m' = 0. This correspondence is valid for both pure [ZOj arid impure 
films |182j . and is given by the formal replacement of the film thickness L$ with 1/T. 

7.6.4- Superfluid-insulator and metal-insulator phase transitions 

The instability of the quantum critical behaviour with respect to the disorder can be in- 
terpreted as a signal for phenomena of localization |2()4[ I2()5|l2()6j . Huang and Meng [204J 
have shown by the methods of the pseudopotential and the Bogoliubov transformation 
that at zero temperature the random impurities can deplete the Bose condensate, though 
not completely. On the other hand, the random impurities generate an amount of normal 
fluid equal to a part of the condensate depletion. 

A paper by Giamarchi and Schulz |2()5j also provides reliable theoretical predictions 
about localization-delocalization phase transitions in disordered superfluids. This paper is 
mainly devoted to the localization problem in one- dimensional electronic systems (metal- 
insulator transition [2*l l207j ) but the analogy with the localization-delocalization problem 
(superfluid-insulator transition |2"H I148j ) in superfluids is also considered. A variant of 
RG is developed |2()5j that is convenient for the treatment of the simultaneous effect of 
interelectron interactions and disorder. The phase diagram of the one- dimensional metal 
is deduced from the stability properties of the RG equations. 
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Remember that the interelectron interaction alone can produce a localization of elec- 
trons and this is the Mott localization. On the other hand, the disorder alone can 
induce the Anderson localization of the eigenstates. The interplay between these two 
localization mechanisms is a matter of an intensive research. The analogy between the 
superconducting-insulator transition in metals and the superfluid-insulator transition in 
Bose superffuids was discussed in Refs. [2U 11481 I^U5j . The superconducting- insulator 
phase transition in disordered thin films and wires was investigated in Refs. |2()8[ 209 . 
As noted in Ref. [21] the outcome of the competition between the interparticle interactions 
and the disorder may result, under certain circumstances, in a gapless insulating "Bose- 
glass" phase. It is supposed that this may happen when the interparticle interaction is 
sufficiently strong [2~4"] . 

A generalization of the Halperin-Lubensky-Ma paper |2()2j for magnetic fluctuation ef- 
fects on the order of the phase transition in pure superconductors was made for the 
zero-temperature superconductor-insulator phase transition in one- and two-dimensional 
pure |21Uj and disordered |211j superconductors (a 5— correlated disorder of type random 
impurities). It has been found that the zero-temperature phase transition is continuous 
and the system behaves like a normal metal up to the phase transition point (the re- 
sistence has a finite nonzero value at T = 0). A recent theoretical study |212j of the 
superconductor-insulator transition in a two-dimensional array of Josephson junctions 
with random couplings demonstrates that the previously predicted Bose-glass phase [2U 
213] is a metal state that has a well defined zero-temperature limit for the conductivity. 

In Ref. |206 one-dimensional Bose system described by a Hamiltonian of type (27) has 
been studied by Monte Carlo techniques with the aim to clarify the competition between 
the strong interaction and the random impurities. The model considered in Ref. [206] is a 
lattice version of the field theory given by Eqs. (27), (28), (89), and (90), only the random 
distribution function is somewhat different. The main result from this Monte Carlo in- 
vestigation is the numerical justification of the Bose-glass phase proposed in Ref. [21] and 
the prediction of a new Anderson-type insulating phase ("Anderson glass") in the weak 
coupling regime. The existence of two insulating phases of glassy type was proposed for 
the first time in Ref. |205j . 

8. Classification of quantum phase transitions 

Here we shall enumerate the main types of quantum phase transitions. The general classi- 
fication in two types: first order (discontinuous) and continuous phase transitions known 
from the classical theory remains valid also for the zero temperature phase transitons. 
Both basic types of phase transitions are possible at zero temperature in accord with the 
Nernst theorem. 

We should have in mind that the quantum effects affect more essentially the close vicinity 
of low- and zero-temperature equilibrium points of first order phase transitions rather 
than low- and zero-temperature critical points of continuous phase transitions. While the 
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correlation length at critical points is infinite and the classical fluctuations dominate at 
any finite (T c > 0) critical point, at the equilibrium points of the first order transitions 
this length is finite and the quantum effects may compete the thermal fluctuations up to 
an infinitesimal vicinity of the transition point. The outcome of this competition between 
the classical and quantum fluctuations depends on nonuniversal properties of the system 
and should be checked for each particular investigation. 

The low- and zero-temperature phase transitions may be either quantum or classical ones, 
depending on whether the quantum correlations have an essential effect on the phase 
transition properties or the classical effects prevail. General thermodynamic and scaling 
arguments (Sec. 2) as well as statistical considerations (Sees. 3 and 6) show that the zero 
temperature limit itself does not guarantee a quantum effect on the zero-temperature 
phase transition. This is valid both for discontinuous and continuous phase transitions. 

Moreover, the zero temperature phase transitions can be classified as T— driven and 
X— driven transitions and the properties of these two types of transitions are quite dif- 
ferent. Note, that there is a wide class of systems as, for example, Bose fluids, which do 
not exhibit X— driven transitions (Sees. 3 and 5). 

The continuous quantum phase transitions, i.e. the quantum critical phenomena can be 
classified in two large groups: 

(1) universal, and 

(2) nonuniversal. 

If all universal characteristic features (critical exponents, including the correction to scal- 
ing exponents) of a quantum system (model) at spatial dimensionality d are identical to 
the respective quantities characterizing the correspondent classical system (model) at spa- 
tial dimensionalities (d+z), where z is the dynamical critical exponent, then the quantum 
critical phenomena in this quantum system are from the group (1). The "correspondent 
classical model" of a quantum field model of the type (27) is that given by setting all 
Matsubara frequencies in Eq. (27) equal to zero. Therefore, the correspondent model is 
the same (the usual 4 -theory) for the whole set of quantum models defined by Eqs. (27), 
(28), (29), and the possible values of the exponents m, and m'. For quantum lattice 
models, one should take the classical limit in the proper way in order to obtain the cor- 
respondent classical model. Within RG, the phenomena of the type (1) are described by 
zero-temperature FPs which are identical in properties to some respective (conjugate) fi- 
nite temperature FP and CQC is the only quantum effect on the zero-temperature critical 
behaviour. 

If the quantum critical phenomena are not in group (1) they belong to the group (2). The 
examples discussed in details in this review are the quantum critical phenomena in Bose 
fluids and XY magnets, described by the Gaussian-like FP (Sees. 5.3-5.4). 



69 



To the best our knowledge this classification scheme is published for the first time in this 
review. The nonuniversality of quantum critical phenomena in Bose fluids and XY systems 
has been pointed out for the first time in Ref. [BE] and then applied in Ref. |lfi3j for the 
interpretation of results for a particular problem of structural phase transitions. With the 
help of various convenient terms of the authors choice the phenomenon of nonuniversality 
is confirmed by researchers in this field as already mentioned in Sec. 5.1. 

9. Concluding remarks 

The following results and notes can be enumerated at the end of our review: 

1. The complete investigation of quantum phase transitions naturally includes the phase 
transitions at low- and extremely low temperatures and cannot be restricted to zero- 
temperature phase transitions only. The quantum phase transtions are a selected part of 
the low-temperature (including zero-temperature) phase transitions. 

2. The properties of T— driven and X— driven phase transitions at low and zero temper- 
ature differ substantially from each other. The quantum effects on the phase transition 
properties are mainly exhibited at X— driven phase transitions at extremely low or zero 
temperature. There is, however, a conformity between the low-temperature T— driven 
continuous phase transitions and the high-temperature X— driven continuous phase tran- 
sitions, as noticed in Ref. [3] and Sec. 2.4. 

3. The first-order quantum phase transitions with sufficiently low equilibrium transition 
temperatures T eq have not been studied enough although they occur more frequently in 
crystal bodies and should be drastically affected by the quantum fluctuations. 

4. The phenomenological scaling theory of continuous quantum phase transitions enters in 
the class of theories of dimensional crossover and takes advantage from similar theories as, 
for example, the scaling theory of finite size systems. Note, that the crossover exponent z/ , 
introduced in Sec. 2.3 can be compared with similar exponents in alternative approaches 
to finite-size scaling and quantum criticality. 

5. Two types of crossover phenomena can be defined in the low-temperature limit: HLTC 
and CQC. The former does not necessarily include quantum effects on the critical be- 
haviour. CQC is the product of quantum effects. 

6. The performance of zero-temperature limit in the scaling equations does not itself 
guarantee quantum effects on the phase transition properties. Rather a phase transition 
can be considered quantum provided the general criterion (3), or, equivalently, the cri- 
terion (9) is satisfied in the transition region. Asymptotic quantum critical behaviour 
occurs unless the respective criterion, (3) or (9), is obeyed in the asymptotic vicinity of 
the critical point. This is quite strong requirement and may be fulfilled in certain sys- 
tems with suitable nonuniversal parameters only at their zero temperature (multi) critical 
points (T c = 0). However, for the finite experimental accuracy, apparent quantum effects 
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may be observed also at sufficiently low-temperature but finite transition points (T c > 0). 

7. The scaling behaviour of basic models of statistical physics does not exhibit quantum 
effects in the zero temperature limit but rather follows the criteria for classical critical 
behaviour. This is valid for certain regimes of performance of the phase transition (Sees. 3 
and 6). 

8. The interacting Bose fluids and XY magnets exhibit a nonuniversal quantum critical 
behaviour. The RG investigations indicate a break down of a well definite form of uni- 
versality in these systems, i.e., the phenomenon of quantum universality breaking takes 
place. 

9. A quantum non-universality occurs also in disordered quantum systems, where it ap- 
pears as an instability of the zero-temperature critical behaviour towards all known mech- 
anisms of disorder. All RG predictions for the critical properties of disordered systems 
at finite temperatures break down at zero temperature where the dimensional crossover 
lowers the critical dimensionality of the fluctuation interactions. 

10. There is a formal analogy between three physically different types of crossover phe- 
nomena: CQC, the (finite-size) crossover from two-to-three dimensional systems, and the 
crossover from point to extended random impurutes in disordered systems. 

11. The available results about quantum phase transitions are not enough to build up a 
satisfactory picture of these intriguing phenomena. 

12. In this review we have tried to show some of the outstanding problems and present 
a unified and consistent framework of the available results. It has been shown that the 
understanding of the quantum phase transitions as "zero-temperature phase transitions" 
is a restricted and incomplete point of view. Besides, the concept of CQC as a main 
result of the quantum effects on the critical behaviour is also unsatisfactory and should 
be complemented by the concept of the quantum universality breaking. 

13. A substantional interest from both theoretical and practical points of view are the 
so-called T— driven phase transitions in the ultra-low temperature scale (T ~ 0). Their 
description is an outstanding problem. Among the unresolved problems are the equation 
of state for quantum phase transitions, in particular, in Bose fluids, where one should rep- 
resent the behaviour of the system in terms of the density of particles and the temperature 
in the quantum limit (T — > 0). 
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Figure captures: 

Fig. 1. (a) Low temperature part of a critical line with zero temperature critical point 
T C (X Q ) = 0. The shaded part (a — X — b) of the transition region, marked by the lines 1 
and 2, corresponds to a low temperature classical behaviour, (b) Low temperature critical 
line with X = 0; domains (a — X — b) and 1 — — 2 coincide. 

Fig. 2. (a) The function f(tp) for a standard second order phase transition for: r = 1.0 
(curve 1), r = 0.5 (2), r = -0.5 (3), r = -1.0 (4). (b) The function /(</?) for the 
IBG at constant density (the curves 1-4 correspond to the same values of r as given for 
Fig.2a). 

Fig. 3. (a) A tree diagram denoting the interaction part of the Bose Hamiltonian. (b) 
The compact self-energy diagram which is equal to zero in the limit T — > [the thick 
loop denotes the full (renormalized) Green function G(q)]. (c) An example of a diagram 
from the perturbation series for the interaction vertex which gives a zero contribution in 
the zero temperature limit, (d) The infinite ladder series of diagrams which yields the 
geometric progression (51). 

Fig. 4. (a) The graphical representation of Eq. (67) for J = 1. (b) The high- and low- 
temperature parts of the curve in Fig.4a with an indication of T— and T— transitions 
(see the text). 

Fig. 5. A scheme representing four relevant FPs (P, R, Gl, Un) in disordered systems. 
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